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A non-Abelian flavour symmetry in a minimal supersymmetric standard model can explain 
the flavour structures in the Yukawa couplings and simultaneously solve the SUSY flavour 
problem. Similarly the SUSY CP problem can be solved if CP is spontaneously broken in 
the flavour sector. In this work, we present an explicit example of these statements with a 
SU(3) flavour symmetry and spontaneous CP violation. In addition, we show that it is still 
possible to find some significant deviation from the SM expectations as far as FCNC and CP 
violation are concerned. We find that large contributions can be expected in lepton flavour 
violating decays, as [i — > ej and t — > /17, electric dipolc moments, d e and d n and kaon CP 
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ON , violating processes as €k- We also show that without further modifications, it is unlikely for 
these models to solve the $s a anomaly at low-moderate tan f3. Thus, these flavoured MSSM 

• ' realizations are phcnomenologically sensitive to the experimental searches in the realm of 

C flavor and CP violation physics. 
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I. INTRODUCTION 

In the next few years, after a long impasse in the phenomenological searches for new physics in 
the energy frontier, all the high-energy particle physics community will be focused on the results 
from the LHC experiments. Its first goal will be the study of the physics of electroweak symmetry 
breaking, but besides this, we also expect some kind of new physics around the electroweak scale. 
Supersymmetry (SUSY) is perhaps the new physics option that is best motivated. The LHC should 
find some SUSY particles if SUSY is indeed the solution to the hierarchy problem and provides a 
candidate for the dark matter observed in the universe. 

Supersymmetry has been extensively studied in the last decades, but most of these studies have 
been done in the framework of the so-called Constrained Minimal Supersymmetric Standard Model 
(CMSSM). The CMSSM is one of the simplest supersymmetric extensions of the Standard Model 
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as it assumes universality of the supersymmetric soft breaking terms and is completely determined 
by four parameters (Mi/ 2 , rno, -Aq, tan/3) plus a sign (sg(/x)). This simplified model is very useful 
to explore the main features of the SUSY spectrum in collider experiments, however, nobody 
really believes that the realization nature has chosen of supersymmetry is exactly the CMSSM, 
specially concerning flavour. In analogy to the known flavour structures in the Yukawa couplings, 
we naturally expect non-trivial flavour structures in the supersymmetry soft-breaking terms. 

Therefore, we should consider other more general flavoured MSSM models as the SUSY models 
we can find when we analyze the experimental results from LHC experiments. Nevertheless, it is 
well-known that the presence of generic flavour structures in the SUSY soft-breaking terms causes 
the so-called "supersymmetric flavour problem". Flavour changing neutral currents (FCNC) and 
flavour-dependent CP violation observables receive too large contributions from loops involving 
SUSY particles and can not satisfy the stringent phenomenological bounds on these processes [l|,Q]- 
Although this statement is still true, it is important to emphasize that the basis of this problem 
lies clearly on our total ignorance about the origin of the observed flavour and CP-violation in our 
theory, and this includes also the SM Yukawa couplings. The real flavour problem is simply our 
inability to understand the complicated structures in the quark and lepton Yukawa couplings, and 
likewise for the soft-breaking flavour structures in the MSSM. It has been recently shown in the 
literature 3|, |4J that an MSSM model with a non-Abelian flavour symmetry allows a simultaneous 
understanding of the flavour structures in the Yukawa couplings and the SUSY soft-breaking terms, 
adequately suppressing FCNC and CP violating phenomena and solving the SUSY flavour problem. 

In this work, we intend to exhibit a concrete example of such strategy considering an MSSM 
(with the limited number of new parameters present in the CMSSM) and showing that the presence 
of a non-Abelian family (horizontal) symmetry can simultaneously account for the solution of both 
flavor problems, namely that it is possible to successfully reproduce the correct fermion spectrum 
while adequately suppressing FCNC and CP violating phenomena. Regarding the phenomenology 
of these flavoured MSSM models, we do not expect large differences from the expectations in the 
CMSSM from the point of view of collider studies, as we know that the departures from flavour 
universality are strongly constrained by the present FCNC and CP violation experiments. Much 
probably the resolution of the new supersymmetric flavour structures will have to rely on FCNC 
and CP violation experiments. In the following we will analyze the FCNC and CP violation 
phenomenology of our flavoured MSSM example. We will show that it is still possible to find 
some significant deviation from the the SM expectations as far as the FCNC and CP violations are 
concerned, making these realizations phenomenologically sensitive to the experimental searches in 
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the realm of flavor physics. 

In the next section, we introduce the flavour symmetries in SUSY and study its effects in the 
SUSY soft-breaking terms. Then, we illustrate these effects with an explicit model with a SU(3) 
flavour symmetry In section [TTTl we analyze the different FCNC and CP violation observables in 
this model including lepton flavour violation, kaon physics, B physics and electric dipole moments. 
Section |IV] is devoted to a combined analysis of all the FCNC and CP violation observables and we 
show the correlations between the most interesting observables in this model. Finally, in section IVl 
we present our conclusions. 



II. FLAVOUR SYMMETRIES IN SUSY 



Flavour symmetries have been used with success in the past to try to extract some meaning 
from the complicated structures of fermion masses and mixings in the Standard Model (SM). Using 
the Froggat-Nielsen mechanism jjj], flavour symmetries explain the structure of the SM Yukawa 
couplings as the result of a spontaneously broken symmetry associated with flavour. In different 
extensions of the SM, these flavour symmetries will also constrain the new couplings and masses. 
For instance, in the context of a supersymmetric theory, a flavour symmetry would apply equally 
to the fermion and scalar sectors. Therefore, this implies that in the limit of exact symmetry 
the soft-breaking scalar masses and the trilinear couplings must be invariant under the flavour 
symmetry and the structures in the soft-breaking terms will be generated after the spontaneous 
breaking of the flavour symmetry. Both the flavour structures in the Yukawa couplings and in the 
soft-breaking terms are generated by the same mechanism and thus we can expect a close relation 
between them 1 . 

The flavour-diagonal scalar masses, i.e. couplings $(f), are clearly invariant under any symmetry 



and are always allowed by the flavour symmetry. However, in general, this does not guarantee that 
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121 ] family symmetry, the 



they are family universal. In the case of an Abelian 
symmetry does not relate different generations, and, generically, different fields will have different 
diagonal soft masses. On the other hand, a non-Abelian family symmetry groups two or three 
generations in a single multiplet with a common mass in the symmetric limit, thus helping to solve, 
in principle, the FCNC problem. This was one of the main motivations for the construction of the 



1 The relation between the soft-breaking terms and the Yukawa matrices is present in gravity mediation mechanisms, 
but it would not be present if the mediation mechanism is flavour blind, as in the case of gauge mediation or anomaly 
mediation 



first £77(2) flavour models Q], where the first two generation sfermions, facing the strongest 
constraints, share a common mass. In the case of an SU(S) flavour symmetry 3, 15, all 
three generations have the same mass in the unbroken family symmetry limit. On the contrary, 
trilinear couplings are completely equivalent to the Yukawa couplings from the point of view of the 
symmetry because they involve exactly the same fields (scalar or fermionic components). Thus, 
they are generated after symmetry-breaking as a function of small vevs. Therefore, to solve the 
"SUSY flavour problem", we will consider in this work non-Abelian family symmetries, and more 
exactly SU (3) theories or discrete versions of this symmetry. 

Apart from these renormalizable mass operators in the Lagrangian, we can construct non- 
renormalizable operators neutral under the flavour symmetry inserting an appropriate number of 
flavon fields. The flavon fields, charged under the symmetry, are responsible for the spontaneous 
symmetry breaking once they acquire a vev. Then, higher dimensional operators involving two SM 
fermions and a Higgs field, with several flavon vevs suppressed by a large mediator mass, generate 
the observed Yukawa couplings. In the same way, these flavon fields will couple to the scalar fields 
in all possible ways allowed by the symmetry and, after spontaneous symmetry breaking, they 
will generate a non-trivial flavour structure in the soft-breaking parameters. Therefore, by being 
generated by insertions of the same flavon vevs, we can expect the structures in the soft-breaking 
matrices and the Yukawa couplings to be related. 

The structure in the Yukawa couplings is not completely determined by the observed values of 
fermionic masses and mixing angles. To solve this problem and fix the Yukawa couplings, we accept 
that the smallness of CKM mixing angles is due to the smallness of the off-diagonal elements in the 
Yukawa matrices with respect to the corresponding diagonal elements, and we make the additional 
simplifying assumption of choosing the matrices to be symmetric. With these two theoretical 
assumptions, and using the ratio of masses at the GUT scale to define the expansion parameters 
in the up and down sector as e = sjmsjm^ and e = y/m c /mt, we can fix the Yukawa textures in 
the quark sector to be: 
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where e ~ 0.15, e ~ 0.05 and the xfj are 0(1) coefficients fixed by the observed values of fermion 
masses and mixings. In the Appendix we show the full structure of the Yukawas, and find the best 



fit values xf 2 ~ 1.7 



> -^13 



0.4, — 1-8, xf 2 — 1-4, x" 3 ~ 2, £23 — 2- In the leptonic sector we will 



follow the same strategy as in Ref. [4] and require unification of charged lepton and down-quark 
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flavour matrices, i.e. we embed our model in a grand unified framework, for instance 50(10), and 
we try to explain simultaneously quark and lepton Yukawas. The detailed structure of the leptonic 
Yukawa matrix is also shown in the Appendix. 

Taking this Yukawa structure as our starting point, we will generate the flavour structure of 
the soft-breaking terms in g eneric non-Abelian SU(3) flavour symmetries or discrete versions, like 



A(27) or A(54) 



18 



19|, 



20 



21 



221 ] . Under these symmetries, the three generations of SM fields, 



both SU (2) ^-doublets and singlets, are triplets 3 and the Higgs fields are singlets. Therefore 
Yukawa couplings and trilinear terms are not allowed by the symmetry. In these models, we add 
several flavon fields. For instance in |^] we have #3, 623 (anti-triplets 3), 63 and 623 (triplets 3), 
while in [itJ we have also #123 and #123 • The symmetry is broken in several steps, first SU(3) 
is broken into SU{2) by the vev of #3 and 03, with (#3) = (0,0, a) and a being of the same 
order as the mediator mass Mf, i.e. a/Mj ~ O(l)- Subsequently 6*23 and 623 get a smaller 
vev oc (0,6,6), with b/M^ ~ e and b/M u ~ e. Notice that in principle, we have three different 
mediator masses, Mf = Ml, M u , M^, because the flavour symmetry must commute with the SM 
symmetry and therefore the vector-like mediator fields must have the SM quantum numbers of the 
usual particles 2 . In some models we also have #123 and #123 getting a lower vev cx (c, c, c) with 



c/M d 



Notice that the mentioned vevs require a vacuum alignment mechanism. In this work we do 



not specify a particular mechanism, but we refer the interested reader to the examples in [3j, 



23j, 



24 



ia 



The basic structure of the Yukawa superpotential (for quarks and leptons) is then given by 



in the absence of #123 fields, or 



e\6{ + {6363) + e M e 2 3,k0zAz (Ma) + • • • 



W Y = 



q + $23^23 4" ^930193 + 0193093 • • • 



'3 U 3 



23^123 



'123^23 
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(3) 



in the models with #123- All flavon fields in these equations should be understood as Oi/Mf. Note, 
however, that the SU(3) symmetry is not enough by itself to determine the required structure 
in the superpotential. Generically, we have to introduce additional global symmetries to forbid 
unwanted terms, like a mixed term 63623, that would spoil the Yukawa structure. Nevertheless, 
the structure in Eqs. ^ and ([3D is quite general for the different SU(3) models we can build, and 



for additional details we refer to [3. Il5l. Il6|. 



2 For simplicity, we take Md and M u = Ml in all our numerical calculations. 



In the same way, the scalar soft masses deviate from exact universality after SU(3) breaking. As 
explained above, (j)\4>i is completely neutral under gauge and global symmetries and gives rise to a 
common contribution for the family triplet. However, after SU(3) breaking, terms with additional 



281 ] . Any invariant combination of flavon 



flavon fields give rise to important corrections 0, q, y , i 

fields can also contribute to the sfermion masses (at least with Planck scale suppression). In this 
case, it is easy to see that the following terms will always contribute to the sfermion mass matrices 
(the presence of the #123 field will depend on the model): 

(Mj)i = ml (51 + ^ \e\.e{ + + 4.4, + + e\ 23l e{ 23 + e 12 ^J 

+ ^(^4 fc 4K^ m %,^23,n)+ ... ), (4) 

where / represents the SU(2) quark and lepton doublets or the up (neutrino) and down (charged- 

lepton) singlets, so that Mf = Ml, M u , From here we see that taking (#123) /Mj, = (#123) /Md = 

e 2 , both models with and without #123 have the same structure up to corrections 0(e 4 ) that will be 

always subdominant in the relevant soft mass matrices in the basis of diagonal Yukawa matrices. 

Notice that in Eq. flU) only the fields that enter the superpotential have their vevs and associated 

□ - 

mediators masses fixed by the e or e parameters. For instance in the discrete model of Ref.[19[, #123 
(the 3-field that in this reference is written as #123) does not enter the superpotential. Therefore, its 
contributions to the soft masses can be suppressed even having a large vev if the associated mediator 
mass is high enough 3 . However, in the continuous version of this model [l7|. the vevs of 6123/Mf 
and 8123/ Mf are constrained to be both 0(e 2 ) by D-flatness and thus they are not dangerous in the 
soft-mass matrices. Moreover, we have to remember that these deviations from universality in the 
soft-mass matrices proportional to flavour symmetry breaking come always through corrections in 
the Kahler potential and these effects will be important only in gravity-mediation SUSY models. 

In the case of the trilinear couplings we have to emphasize that from the point of view of the 
flavour symmetry these couplings are completely equivalent to the corresponding Yukawa coupling. 
This means that they necessarily involve the same combination of flavon vevs, although order 
one coefficients are generically different because they require at least an additional coupling to a 
field mediating SUSY breaking (in general coupled in different ways in the various contributions). 
Therefore, from our point of view, we expect that the trilinear couplings have the same structure 
as the Yukawa matrices in the flavour basis. However in general they are not proportional to the 
Yukawas, because of different O(l) coefficients in the different elements. Thus, we can expect that 



3 We thank G. G. Ross and I. de Medeiros Varzielas for clarifying this point. 
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going to the SCKM basis does not diagonalize the trilinear matrices. In fact, the trilinear matrices 
maintain the same structure as in the flavour basis and only the O(l) coefficients are modified. 

Therefore, we can see that in these SU (3)-like models with the three generations unified in a 
single field, we have basically the same "leading order" structures in the soft mass matrices and the 
trilinear couplings directly related to the structures in the Yukawa couplings. In the following we 
concentrate in the SU(3) model of because it has a complete phase structure in the flavon vevs 
consistent with the CKM phase. Remember however, that the main features of the soft terms are 
similar in other models and in principle the phase structure could be also adapted in these models. 



A. SU(3) Flavour Model 



Let us now specify more explicitely the SU (3) flavour model that we use as our main example. 
The full superpotential is determined by SU(3) and several global symmetries are used to forbid 
unwanted terms that would spoil the observed structure of the Yukawa couplings. Using the charges 
presented in Table HU the leading terms in the superpotential are, 

6\6{ + #23#23 S + ( e * fc ^23,fc^3,Z#23 + ^^23 ^3,1® 23) (#23^3) 



Wy = H^tf 



+ 



e J,?7 #23,z (^23^3) + ^ l &3,l (6*23^3) (#23^23) + • • • 



(5) 



where to simplify the notation, the flavon fields have been normalized to the corresponding mediator 
mass, which means that all the flavon fields in this equation should be understood as 6i/Mf. The 
field £ is a Georgi-Jarlskog field that gets a vev in the B — L direction, distinguishing leptons and 
quarks. Furthermore, as said above this model is embeded in a 50(10) grand unified structure 
at high scales, which allow us to relate quark and lepton (including neutrino) Yukawa couplings. 
However, the SU(2)r subgroup of 50(10) must be broken as we need different mediator masses 
for the up and down sector and, in fact #3 and #3 are 3 © 1 representations of SU(2)r which is 



broken by their vevs 
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After spontaneous symmetry breaking the flavon fields get the following vevs: 

M , x M 



(ft 





V 1 / 

hz) 



of e ix 







a? e 



d ia d 



a, e 







\ 



&23 

V ^23 e ift J 



?23/ 



b 23 



\b 2 3 



(6) 



8 



Flavon Phase 




X 


Pa 




Allowed Values 


(-40 ±2)° 


(20 ± 10)° + 180°n 


-(20 ± 10)° + 180°n 


Unconstrained 



TABLE I: Flavon phases after imposing CKM constraints, with n E N. Notice that the approximate errors 
we quote in \ and P3 are highly correlated, as shown in the Appendix. 



where we require the following relations: 
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(7) 



These relations are valid at the flavour breaking scale, that we take as the GUT scale in the 
numerical evaluation. 

It is straight-forward to see that this superpotential reproduces correctly the required Yukawa 
structure in Eq. ([T]). For completeness, we also list in Table U the values that the flavon phases can 
take, given the constraints imposed by the CKM matrix. The analysis leading to such constraints 
is found in the Appendix. 

We can now turn to the soft breaking terms. As mentioned in the previous section, a universal, 
flavour diagonal mass term will always be allowed. Moreover, in a SUSY theory, the same messenger 
fields as in the Yukawas will couple the flavons to the scalar fields in the soft terms. Thus, the 
e and £ parameters still act as expansion parameters, and represent important corrections to the 
soft terms. 

Clearly any coupling involving a flavon field and its hermitian conjugate (i.e. is invariant 

under the flavour symmetry. From this we can deduce that the soft mass terms get a minimum 
structure determined uniquely by the flavon content of the model and on their vevs. This minimum 
structure is obtained from the following effective terms: 



(Mf)i = m3K4 



OlM + 03,i03 + 023j023 + 023,i023 



7T Tit j 



+ (^0 3 %)(e jmre 0L4j + °23 )(^03, m 023,n) + 



t fc 7Tt 



(8) 



In Table |TT] we show a choice of global charges that reproduces the correct Yukawa structure and 
does not allow other terms at leading order in the Kahler potential (soft-masses). 

In the squark sector, after rephasing the fields such that the CKM matrix elements V uc i, V us , 
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TABLE II: Charges required to build the minimal RVV1 Model. 



V c b and Vth are real, the soft masses in the SCKM basis are: 
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where M? (M~) is in the basis where (Y e ) is diagonal, uj us is related to the CKM phase, and 
lo' = uj us — (5 U — Sd)- The uj us and Si phases can be found in the Appendix. The structure of M~ 
in the basis where Y u is diagonal is similar to M 2 C ^. Notice that, although the structure in terms 
of e, e in M~c and M~ is the same as that of M~ and M? respectively, the coefficients from the 

e R L d c R Q 

SCKM rotation and RGE evolution are different due to the Georgi-Jarlskog field X. This can be 
seen in the Appendix, and is summarized in Table [V] in Section IIII Al The phase structure of the 
slepton matrices is different to the one of squarks, since the latter have been rephased in order for 
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TABLE III: Charges for RVV2 Model. 



the CKM matrix to follow the standard phase convention. 4 

In Eq. ([9]), we have written only the leading contribution in e and e to each element with a 
leading phase, omitting effective complex 0(1) coefficients. This sets the size of the modulus of the 
mass-insertion. The 0(1) coefficients are the remaining part of the full element, after factorizing 
the terms explicitly written in Eq. (|9|). There are many contributions to these coefficients: first, 
we have the real O(l) constants in front of each term in the Kahler potential, at Mqut- We also 
have contributions from the real O(l) constants in the Yukawa matrices, coming from the rotation 
into the SCKM basis. A third contribution comes from the RGE evolution to the EW scale. And 
finally, it is possible to have further contributions from subleading terms in the Kahler potential, 
still at Mqut- 

All of these contributions can involve the flavon phases, so the effective 0(1) can be complex. As 
we are factorizing the leading phases, the phase in each coefficient shall only appear in subleading 
terms. If the leading phases cancel for a particular observable, these subleading phases shall be 
important. Such a cancellation can happen in EDMs. For example, in mass-insertion notation 



30, 



311 ] . the corresponding effective 0(1) coefficients of (Sf s )jm and (5f 3 )LL have the following structure: 

2t(x-A0 



)RR 

Qfs L 



1 - Vb e" 
l-( £ 2 /e 2 ) e - 2 «3). 



(10a) 
(10b) 



Although Eq. (jSJ) is the minimal structure present for all possible models, it is possible, for 
particular choices of the global symmetries and charges, to build other symmetry-dependent soft- 
mass structures. In fact, the observed structure in the Yukawa couplings does not fix completely 
the introduced global charges and it is possible to add new invariant combinations of flavon fields 
to the Kahler potential without modifying the Yukawas. 

The first example of these new combinations of flavon fields in the Kahler is achieved by allowing 



4 As physical observables are independent of phase conventions, we could also rephase the slepton superfields and 
use the same phase structure as squarks. 
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a ^3^23 term (RVV2). The required charges are shown in Table llXH It is easy to check that, with 
these charges, the structure of the Yukawa couplings in the superpotential remains unchanged. 
This is due to the fact that the superpotential is a holomorphic function of the fields while the 
Kahler is only a real function. When rotated to the SCKM basis, the soft-mass matrices become: 
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One can see that the effect of this term in m^ c is to exchange one power of e by a yjj suppression 
in (Sf 3 )nR and (S^rr- In m 2 - , the same terms change an e 2 by an e y ' 5 . However, for tan /? = 10, 
and considering that ewe 2 , such replacements leave the structure of the mass matrices very similar 
numerically to the original one (notice that yt and y^ are taken at Mqut)- Nonetheless, it must 
be remarked that the phase structure of the whole mass matrix is modified. 

As in RVV1, it is crucial to take into account that relative phases exist within the effective 0(1) 
coefficients. For instance, although the global phase of {5\2)aa is still w ns , the 0(1) structure of 
(5f 2 ) LL in RVV2 is now: 



LL\RVV2 



(Ql2 ) 



1 



70 (i + 



J(2 X -0 3 -P' 2 ) 



)■ 



(12) 



Notice that the factor (e 2 /e) =0.45 does not really provide any suppression at all. This means 
that the imaginary part of (5^ 2 )ll is larger than just e 2 esmu us . We label this new, larger, effective 
phase as u>' us . 
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Field 


v 


-0 C 


H 


£ 


o 3 


#23 




#23 


SU(3) 


3 


3 


1 


1 


3 


3 


3 


3 


U(l) 


-1 


-1 





5 


1 


-2 





6 


U'(l) 











-1 








1 


-2 



TABLE IV: Charges for RVV3 Model. 



A second possibility is to allow a (^ lkl ^9 l 23 ^ 9% term (RVV3) in the Kahler, with the charges 
being shown in Table HVl The soft matrices, when rotated into the SCKM basis, have the following 
structure: 

/ 1 _i_ ^2 , 



Mi 



T 



Mi 



M 



Q 



Mi 



Mi 



\ eytt 



l + e z Vt 
-e 3 e~ iw ' 1 + e 2 



e 3 e iJ E y t _ e i(uj'-2x+l3 3 -6 u ) \ 



l + e 2 y b 



_g3 g— iw ua 



\zybe 



-i(uua+/33—5d) 



-e 3 e luJus e y b e 
1 + e 2 
e 2 



e 2 e~ 21 * 

i + yt 

(u>us+l33-S d ) \ 
2 



17ln 



[13a) 



£ 

1+2/6 



run 



(13b) 



-ee yt.e 



i(u) ua -2x+P3+5 d ) 



zyte 



-££ 2 y t e~ i ( UJua ~ 2x+ P 3+Sd " > 



\ zyte 



-i(u) us -2x+/33+S d ) 



( 



l + e 2 y b 

?3 



-e° 
1 + e 2 



1+e 2 
e 2 yt 



i(x-s d ) \ 



e 2 yt 



i + y t 



ml 



J 



(13c) 
(13d) 



/ 



V 



-ee 2 y t e^ 2 *^-^) e y t e -i(x-*d) \ 



-ee 2 y t e *(2x-/%-« 

J(x~Sd) 



E 2 yte -i(x-p3) 



e 2 y t e i( - x ~M 



i + yt J 



(13e) 



ey b e l 

g 2 e i(x-ft) 

\ey b e~ i ^- 5 ^ ^ e -i( x -/3 3 ) j + yb j 
l + e 2 y t 

1 + e 2 

eyte 

with 5j defined in the Appendix. 

This model shows larger deviations from RVV1 in the LL sector. It is important to notice that 
(^12) 'll is now of order ee 2 instead of e 2 e, which will have considerable consequences in processes 
such as fx — > ej. Likewise, (df^LL is of order eyt instead of e 3 , so an enhancement in r — ► e-y 
should be expected. Regarding the RR sector, for tan/3 = 10, the y b suppression at Mqut has 
roughly the same size as an e suppression, so once again the structure of m 2 ~ is numerically similar 

d R 

to RVV1 for this value of tan f3. 
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The trilinear couplings, on the other hand, follow the same symmetries as the Yukawas. Thus, 
they have the same flavon structure in RVV1, RVV2 and RVV3. Nonetheless, although they have 
the same structure, they do not have the same 0(1) constants, which means that the rotation into 
the SCKM basis does not diagonalize them. In the SCKM basis, after rephasing the fields, the 
trilinears have the following structure: 

( Q £ 3 e -iw' £ 3 e i(2 X -u') \ 

e 3 e iuJ ' e 2 S u e 2 S u e 2i * 

( 



A u 



Ao yt 



(14a) 



e 3 e _<IJu 



e 3 e~ iUu 



gr3 g — ibJu 



e 2 E d e 2 T, d 



An 



(14b) 



I ^3 e i(a; us +2/3 3 -2x) ^2 ^ e 2i((3 3 - x ) i j 

where we have neglected the O(l) coefficients. One can get A e by taking A d with (E e ) = 3 (Ed). 

The soft mass matrices in Eqs. (I9|ll|13fl"4"l) are given at the large scale ~ Mqut- Then, we 
have to include also effects coming from the running from Mqut to M\y. Although these effects 
produce further non- universal contributions, they are usually smaller than the terms presented 
here. In the quark sector, the misalignment of the Y u and Yd matrices gives sizeable contributions 
to the LL and LR sectors, analogous to the MFV contributions of CMSSM models. In the lepton 
sector with RH neutrinos, the same happens due to the misalignment of Y v and Y e 32 



33 



3J. 



Both of these contributions are unavoidable, albeit the Y v contribution is highly model dependent. 

Moreover, in the present case there are new contributions to the running given by the intrinsic 
non-universality of the soft mass matrices. Although these effects contribute even in the RR sector, 
it turns out that the magnitude of the generated off-diagonal terms are, at most, of the same order 
as those at Mqut- This means that the addition of running effects will only change the already 
unknown 0(1) constants, such that the low-energy phenomenology can still be understood by 
analyzing Eqs. (^|). (fTTj). f[T5|) and ([13]) . 

Finally, we have to emphasize once more that all the soft matrices presented here have unknown 
O(l) coefficients in all non-universal terms. The flavour symmetry allows us to fix the order in e 
and £ of the different entries but the final values could easily vary up or down by factors of two 
and this has to be taken into account when analyzing our numerical results. 
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III. FCNC AND CP VIOLATION OBSERVABLES 

In the following, we shall study flavour and CP violating observables, presenting analytical 
estimates in mass insertion approximation (MIA). Moreover, we will perform a full numerical 
analysis in the SUGRA parameter space through a scan in uiq and Mi /% for fixed values of tan/3 
and ao = A^/m®. The numerical analysis is done defining the Yukawa, trilinear and soft mass 
matrices at Mji av = 2 x 10 16 GeV, for the three different versions of the model, as explained in 
Section III Al Then the different flavour matrices are evolved to the electro- weak scale, solving 1- 
loop RGEs with SPheno [sslj . O(l) coefficients in the Yukawa matrices are determined by requiring 



361 ] . The result of such fit is presented 



a good fit on the fermion masses and quark mixings at Mz 
in the Appendix. 

Regarding the unknown O(l) constants in both the superpotential and Kahler potential, as we 
do not have a full high-scale model, we shall fix each 0(l)s at a random value, between 0.5 and 2. 
Notice they can be of either sign. Thus, a model is characterized by both the choice of symmetries 
involved, and by the particular 0(l)s we have in front of each effective term. In Section ITVl we 
shall take into account their variation. 

After running the resulting matrices down to the Mz scale, we diagonalize the Yukawas in order 
to obtain the left and right mixing matrices and rotate the soft matrices into the SCKM basis. At 
the Mz scale, for each point of the SUSY parameter space, we compute the SUSY spectrum and 
check that the electroweak symmetry breaking does take place and no tachyonic particles arise. 
Moreover, to be conservative, we require that the Lightest Supersymmetric Particle (LSP) is the 
lightest neutralino. 

In this work we do not include the non-holomorphic corrections to the Yukawa couplings, so 
our results are valid for low and moderate values of tan j3. In this regime, the effects of these non- 
holomorphic corrections are usually small. The only exception is found in the case for the neutron 
EDM, where new imaginary parts induced by these corrections can give sizeable contributions. We 
give details on how we introduce such corrections to this observable in Section IIIID1 

Finally, we apply the following constraints: 



Bounds on the sparticle masses from direct searches at LEP and Tevatron [371 ] . 



Lower bound on Higgs masses from LEP 



. Although the SM bound on the Higgs mass 



is of 114 GeV, in SUSY the modification of the lightest Higgs coupling could reduce this 
constraint. In order to take this effect into account, for each point we calculate the Higgs 
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mixing angle a and use Eq. (17) of Ref. 39J to estimate the correct Higgs mass bound. We 



then use SPheno to calcu 
uncertainty of 3 GeV [40 . 



ate the two loop Higgs mass, and take into account a theoretical 



4l(. 



b — > sj. The experimental world average from the CLEO 



421, Belle 



collaborations is given by 



43[ and BaBar [44] 



BR(B X s 7)^>i. 6G eV = (3.52 ± 0.23 ± 0.9) x 10" 4 . (15) 
We have to compare these values with the MSSM predictions. In our numerical calculation, 



we use the expression presented in Ref. |46l ] in which the branching ratio is explicitly given in 
terms of arbitrary complex Wilson coefficients CV and Cg- At present, the SM contribution 
to this decay is already available at NNLO while the SUSY contribution in a general MSSM 
is partially known at NLO. In this work, we include the NLO SM contribution with a 
modified low value of the scale for the charm mass to reproduce the NNLO SM contribution 
(BR(5 -» X s -f) SM = 3.15 x 10~ 4 ) 



47 



IS 



49J. We add the super symmetric contributions at 



one loop, and require that the total result does not deviate from the experimental value of 
Eq. (|15j) in more than 2-sigma. 



Muon anomalos magnetic moment, a^ = (g — 2)^/2. At present, the experimental result for 
this observable is given by [50I ]. 



oxp 



11659 2080(63) x 10 



-11 



(16) 



while, computing the hadronic contribution by means of the hadronic e + e annihilation 
data, the SM theoretical expectation is 



51 



52, 



53j, 



a s ™ = 11659181(8) x 10 



-10 



The resulting discrepancy is: 

Aa, 



a exp _ a SM 



+302(88) x 10 



-11 



(17) 



(18) 



It is well-known that in the MSSM receives contributions from an d loops 54]. 

Such contributions are approximately given by the following expression: 



a^ SSM / tan p \ ( 300 GeV \ 2 / \iM 2 

1.5 



(19) 



1 x 10~ 9 V 10 J V rn„ ') \ml 

A comparison with Eq. (| 18j) implies that the present discrepancy strongly favours the \i > 
region of the SUSY parameter space. In the case of the theoretical prediction based on r 
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mM m L )is\ \(S e LL U 


KSrvU \(s e RR U \(5 RR U 


RVV1 
RVV2 
RVV3 


\e 2 e y t e 3 3y t e 2 
3y t ££ 2 y t £ 32/ t e 2 


1=3 1=3 =2 

gO gO t 

±e 3 y b e e 2 



TABLE V: Order of magnitude of LFV mass-insertions, for the three models. 





Present Bound 


Future Sensitivity 


BR(/j, -> ei) 
BR{t -► ni) 
BR(t -> e7) 


1.2 x 10" 11 [56] 
1.6 x 1CT 8 [58, 59] 
1.1 x 10~ 7 [61] 


O(10~ 13 ) [57] 
0(l(r 9 ) [60] 
O(10" 8 ) [62] 



TABLE VI: Present bounds and future experimental sensitivities of lepton flavour violating processes. 



decay 



551 ]. the difference of Eq. (|18f) is reduced to ~ la 



521 ] but it still requires a positive 



correction and disfavours strongly a sizable negative contribution. 



A. Lepton Flavour Violation 

As pointed out in jj], flavour models based on SU(3) give rise to potentially large rates of LFV 
processes, such that positive signals of LFV can be found in the currently running or near-future 
experiments, at least for SUSY masses within the reach of the LHC. The arising of large mixing 
among flavours relies on the features of the SU(3) model discussed in the previous sections: the 
presence of nonuniversal scalar masses already at the scale where the SUSY breaking terms appear, 
and the fact that the trilinear Af matrices are in general not aligned with the corresponding Yukawa 
matrices. Let's start considering the case = 0, where the latter effect is strongly reduced so that, 
in terms of mass insertions, BR(/j — ► ^7) mainly depends on \(S1 L )ij\ 2 and | (<5fj/j)r/| 2 . Looking at 
the structure of the slepton soft mass matrices in the three versions of the model (Table |V|) , we see 
that RVV1 and RVV2 are expected to give similar predictions for BR(// — > cy) and BR(t — ► ^7), 
with possibly sizeable contributions coming from both the LL and the RR sector. In the case of 
RVV3, the prediction for BR(r — > ^7) will be also similar to the previous two cases, while we 
expect BR(/i — ► ej) to be strongly enhanced. In fact, for RVV3, the LL mass insertion is larger by 
a factor 9yte/e = O(10) with respect to RVV1 and RVV2, and the BR(/j — ► cy) is consequently 
increased by two orders of magnitude. 

To summarize, let's compare the expections for the different LFV processes. In the case Aq = 0, 
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considering for simplicity only the contribution from 8% L , we have: 



0(1) (RVV1, RVV2, RVV3) (20) 



BR(n^e-y) \m^J T T (5£ L )? 2 

O(10 3 ) (RVV1, RVV2), O(10) (RVV3) (21) 



BR{n^ ei ) ~ V"W r r (5£ L )f 2 

where (r T ) is the // (r) full width. Given the present limits and future sensitivities of LFV 
processes shown in Tab. IVIl we see that BR(r — > ej) is not able to constrain the parameter space 
better than BR(^ — > e~f) in none of the three models. On the other hand, we expect from Eq. (|2ip 
that the present constraints given by \i — > e'y and r — > fi'y, that differ by three orders of magnitude, 
are comparable for RVV1 and RVV2, while fi — > cy should give the strongest constraint in the case 
of RVV3. 

In the case Aq ^ 0, generally large 5 e LR insertions arise as a consequence of the misalignment 
between Aj and the corresponding Yukawa matrix Yf. In this case, the neutralino contribution to 
BR(/x — > ej) gets strongly enhanced J] and the present (or future) bound requires heavier SUSY 
masses to be fulfilled, specially in the region where the gaugino mass is much larger than the 
common sfermion mass. Nevertheless, we expect this effect to be visible only in the case of RVV1 
and RVV2, while for RVV3 the very large (5|,x)i2 should still give the dominant contribution. 

Let us now consider the results of the numerical analysis for the LFV decays. After fixing the 
unknown O(l) parameters to random values, we present in Fig. [I] the current bounds provided by 
fi — > e7 and r — ► /X7 in the (mo, -M1/2) plane, and also the final reach of the MEG experiment. In the 
first row, the Aq = case is displayed (for tan/3 = 10). We see that, as expected, RVV1 and RVV2 
give similar results both for BR(^ — > cy) and BR(r — > fi'y). The regions of the parameter space 
excluded by the present bounds BR(/x — ► &y) < 1.2-10 -1 , BR(r — ► ^7) < 1.6-10 -8 are comparable, 
as obtained by the naive estimate in Eq. (|2ip . although r — > fi'y turns out to be more constraining 
for mo > Mi/2- The presently allowed region is approximately (mo, Mi/ 2 ) ^ (700, 300) GeV. In 
the case of RVV3, fi — > cy already gives a strong constraint, (mo, M 1 / 2 ) > (1400, 800) GeV, which 
is much more stringent than the one provided by r — > fi'y. Although the precise values may vary 
with different O(l) parameters, these plots give an idea of the reach of these LFV experiments. 
As a consequence, for SUSY masses lying within the LHC reach, RVV3 results already rather 
disfavoured, while RVV1 and RVV2 are not strongly constrained. Considering the sensitivity 
expected at the MEG experiment for BR(^i — > cy), O(10~ 13 ), we see that also RVV1 and RVV2 
will be tested in most of the parameter space accessible to the LHC, while RVV3 will be completely 
probed well beyond the LHC reach. Moreover, in case of larger values of tan/3 (e.g. tan/3 = 30), 
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FIG. 1: Current LFV constraints in the m -M 1/2 plane, for RVV1 (left), RVV2 (center) and RVV3 (right) 
for tan/3 =10 and Aq = (first row), Aq = mo (second row). The brown region corresponds to the present 
limit BR(/i — > cy) < 1-2 • 10 -11 , the light brown region represents the sensitivity of the MEG experiment 
(1CP 13 ). The present BaBar+Bcllc combined limit BR(r — > fi-f) < 1.6 • 10~ 8 is shown in dark green, and the 
light green region corresponds to the sensitivity of a Super Flavour Factory (10 -9 ). The dark brown region 
show areas excluded by having a charged LSP or by LEP, excepting the Higgs mass bound, which is shown 
in thick dashed red lines. 



since BR(/i — > e^f) oc tan 2 (5, MEG will be able to test all the displayed parameter space also for 
RVV1 and RVV2 (cfr. [4], where a similar model has been studied). 

In the second row of Fig. [H the three versions of the model are displayed for the case Aq = mo, 
in order to show the potentially large flavour mixing induced by the A-terms. As expected, we 
see that, for RVV1 and RVV2, the [i — > ej bound is increased with respect to the ^0 = case, 
specially for M 1 / 2 ^ rno, as a consequence of the large (<5£/j)i2 insertion contributing in diagrams 
with pure B exchange. Also in this case, MEG has a very high capability of testing the parameter 
space. Indeed, for moderate slepton masses the neutralino contribution is so large that the models 
will be fully probed. Only in the case of rather heavy sleptons, the BR(/i — > e^f) can be suppressed 
enough to escape the reach of MEG. In the case of RVV3, the contribution from LL insertion 
remains dominant and no substantial changes are observed with respect to the case of vanishing 
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trilinear terms. 

As RVV3 is heavily constrained by LFV, in the following we shall exclude it from our analysis, 
and concentrate exclusively on RVV1 and RVV2. 



B. Kaon physics 



Supersymmetric contributions to K°-K mixing are mainly determined by the mixing between 
the first two generations of down type squarks. In terms of mass insertions, a major hadronic 
constraint on ^sm{(Sf 2 )AB} comes from ex- The kaon mass difference, AMk, does constrain 
| (612) AB |- However, this observable is less sensitive to small Mis since the long-distance effects are 
difficult to estimate. 

On the experimental side, current data [3?! shows that |e^ p | = (2.229 ± 0.012) x 10 -3 and the 
phase of €k is <fi e = (43.51 ± 0.05)°. On the theoretical side, many improvements are being made 
in the determination of input parameters needed to calculate the SM prediction of ex- We have 
good information on the CKM matrix elements thanks to the vast amount of data from B physics. 
Also available are new lattice estimates of the Bx parameter which enters the AS" = 2 hadronic 
matrix element. However, these developments have lead to a puzzle in understanding the data. 
The authors of Refs. 



63, 



641 ] realized that within the SM, the theoretical prediction of ex might 



not be enough to account for the measured value if one accepts the unitarity triangle (UT) fit using 



B physics observables. One can summarize their result in the numerical form 

2 



63. 



\ e K I 



1.78 • lQ-*x 



0.92, 



B K \ f \V cb \ 
0.72 J V0.0412 

2 



Rt 
0.914 



0.74 x 



\Vcb\ V ( Rt \ / sin 2(3 
0.0412 J VO-914 J V 0.675 



sin/3 
0.362 



0.40 x 



Vet \ ( S Q (x c ,x t ) 



0.47 



0.00221 



Vtt \ ( Sp(x t ) 
0.5765/ V 2.30 

Vcc 



+ 



0.14 x 



1.43/ V0.000250 



(22) 



The definition of each symbol and its value can be found in the above references. Notice that the 
sin (3, sin 2f3 appearing in this equation refer to the UT angle (3, not to the ratio of the Higgs vevs 
in SUSY. The denominator of each fraction is the central value of the input parameter appearing 



631 ] . Major uncertainties are those in Bk, \V c b\ 4 , Rf, which are 



in the numerator, quoted in Ref. 
5%, 11%, 7%, respectively. One important refinement made in this expression, which had been 
overlooked in most of the literature, is the factor n e that parameterizes suppression of the result 
due to <fi e / 7r/4 and the imaginary part of the 0-isospin amplitude of K — > -kit. In addition, latest 
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attice values of Bk are lower than previous determinations. Here we take Bk = 0.720(13)(37) 



65, 



66 ] (see also 67j]). These two facts cooperate to make |ejj | insufficient for the central values, 
obtaining |eff*| = 1.78 x 10~ 3 , which is smaller than the observed \zk\ by 18%. This deficit should 
be compensated by some new source; since our model has extra sources of CP violation, we are 
lead to invoke these new sources to resolve the puzzle. 

Before we proceed further, a clarification is in order regarding the ex puzzle. It stems from the 
tension among different observables used to determine the UT on the (p,rj) plane. The regions 
preferred by ex, S^k, and AMs a / AMs d do not precisely overlap with one another (see Ref. [68 1 
for example). As we said, if one accepts the regions determined by S^k and AMbJAM^, then 
€k should be modified. Conversely, one may also conceive a scenario where the ex region is in 
fact correct whereas one or more of the B sector observables have been contaminated by new 
physics effects. This could be an equally legitimate solution in a general case. However, it is not 
viable within our framework. The reason is that supersymmetric contributions to the above three 
observables are (up to O(l) uncertainties) correlated and that the ex region is modified much more 
than the other two. We shall add more quantitative comments on this in the next subsection. 

In order to evaluate the supersymmetric effects on e^-, we should calculate M|q| y = 
(K\H^q SY \K). For this, we follow Ref. 69( which presents the AS = 2 Wilson coefficients from 
the gluino-squark box graphs (which are the largest SUSY contributions in the presence of sizeable 
squark mass insertions) and the expressions for evolving those Wilson coefficients from the sparticle 
mass scale down to the hadronic scale. After that, the SUSY contribution to ex is given by the 
formula, 

Im(M|Yf) 



^SUSY 
e K 



sin, 



AM K 



(23) 



And then we can express the supersymmetric contribution to ex in terms of mass insertions as, 



C SUSY 
e K 
exp 
e K 



V2 



sin 



Im< 



(S 



d \2 
12 ILL 



+ (6[ 



d \2 
12 1 RR 



B K 



0.0062 2 



0.72 



2.2 



0.00013 2 



0.70 



0.05 



k{x) \ 
(-1/30) ) 




where x = m|/m~ and functions fe(x) and fe(x) can be found in Ref. 
term proportional to B 5 (fj.) which happens to contribute only 0(1%). 

In the numerical analysis, we set Br = 0.72 
from Ref. Yh 



(24) 



69(. We have omitted the 



65, 



661 ] and use the central values of £?2,...,5(/-0 
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FIG. 2: Contours of |e|^ SY | for RVV1, in the m -Mi/ 2 plane for tan/3 = 10 and A = 0. The lines are (from 
thicker to dotted) |e|^ ISY | = 10~ 3 ,5 x 10 _4 ,10 _4 ,5 x 10~ 5 , respectively. We show predictions for positive 
(left) and negative (right) 0(l)'s for M| c , giving in some cases ef¥ with the wrong phase, which are shown 
in grey. Current LFV bounds are also shown in green, and the (g — 2)^ favoured region is shown hatched 
in yellow. The dark brown region show areas excluded by having a charged LSP or by LEP, excepting the 
Higgs mass bound, which is shown in thick dashed red lines. 

In our flavour models, RVV1 and RVV2, we have \(df 2 ) RR \ ~ ^ - °- 003 and I0*uj)ll| ~ 
e 2 e ~ 0.0004 with different phases in the two models, uj us ~ 0(e) and uj' us ~ O(l) respectively, 
as can be seen from Eq. ([9]) and the paragraph following Eq. (|12|) . Therefore, we can expect that 
the supersymmetric contribution to €k can easily be comparable to e^ p , in particular from the 
contribution in the second term in Eq. (|24p . The additional contribution might well be just what 
we need to fill the gap between and e^ p , which could amount to 18% or more. In Figs. [21 we 
present contours of ex on the (mo,M]y 2 ) plane for RVV1. The left and the right figures are for 
positive and negative 0(l)'s for M~ c , respectively. 

The phase of e|^ SY can be either 43.51° or —136.49° which are respectively marked in black 
and gray in Figs. [2j Note that the phase of is 43.51° and that e|^ SY and should interfere 
constructively to fit the experimental value. Therefore, gray contours worsen compatibility between 
the theoretical and the experimental values of ex by increasing their discrepancy. The phase of 
£ SUSY - g determined by the sign of Eq. (I24p . It is often dominated by the term proportional to 
lm[(df 2 ) Rr\- I n this case, the phase of e^ SY can be shifted by ir if one multiplies either 
{$12) ll or (^12) RR by —1 (of course, the modulus |ef^ SY | does not remain exactly the same due 
to the first term). For example, we can flip the sign of the RR insertion by changing the signs of 
the 0(1) coefficients in M? . This is the reason why the region with gray contours on the left plot 

d R 

turns to black on the right. 

For RVV2, the phase in (Sf 2 )LL is much larger than that in RVV1, as we can see in Eq. (|12p . This 
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makes the SUSY (<5f 2 ) rr(£ '12) ll contribution to e|^ SY dominant in all of the evaluated parameter 
space, so only one of the two possible phases for (5f 2 )iiR solves the ex tension. Furthermore, each 
of the contours in Figure [2] would be scaled roughly by a factor 1/e ~ 6. 

The requirement that the supersymmetric contribution solves the ex puzzle would define a strip 
in this plane which will be shown in all of the following plots, for both RVV1 and RVV2. One 
virtue of such a strip is that the parameter space therein will lead to a definite correlation between 
observables which are not apparently related to each other. We will elaborate on this in Sec. IIVI 

Another important quantity in kaonic CP violation is e'/e. It is specially sensitive to chirality- 
flipping mass insertions, (Sf 2 )LR and (Sf 2 )RL 3QJ]- Using the formulae in Ref. 7JJ], one can write 



the supersymmetric contribution to e'/e in the form, 



/,/\SUSY 

^ = Im 



($12) lr - (5 12 )rl 
1.6 x 10" 5 



xB G 



500 GeV \ 
m- q J 



(25) 



which has been derived for m~/m~ = 1. The 0(1) constant B G parameterizes the uncertainty in 
the hadronic matrix element of the chromomagnetic operator. We have suppressed other factors 
that can change by 0(1) for m~/m~ 7^ 1, since the following argument depends only on the order 
of magnitude of e'/e. 

In this model, there can be flavor-violating A-terms arising from mismatch of the O(l) co- 
efficients between the Yukawas and the A-terms. To estimate their possible effects, we define 
Aq, = clq ttiq, with Aq being the overall dimensionful coefficient of the A-terms. Using Eq. (|14b ) 
and (|9b), one can get rough relationships among different 1-2 mass insertions, 

(Sf 2 )LR ~ (Sf 2 )RL ~ a ^ x (Sf 2 ) RR ~ a x (Sf 2 ) LL , (26) 

where nib is the 6-quark mass at the GUT scale. If we require that the size of the second term in 
Eq. (|24p does not exceed unity, we find that the extra contribution to e'/e has an upper limit like 

jy^SUSY 



< 8.1 x - a ^, (27) 
e mo 



(g')exp 

for purely imaginary (Sf 2 )iR or (Sf 2 )RL- in this case, the supersymmetric fraction within e'/e 
should be lower than 40% for a reasonably high m,Q > 100 GeV. In fact, the phases of these flavor 
and chirality changing insertions are ±w us which are smaller than maximal, as can be seen from 
Eq. (|A5p . Thus their effects on e'/e are smaller than the above estimate. We have numerically 
checked that the gluino-squark loop contribution to e'/e is less than 15% on the strip compatible 
with ex. Given the large theoretical uncertainty in e'/e, this amount of contamination by new 
physics should be hard to disentangle. 
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Another process that could, in principle, be interesting is the decay Kl — ► k^vv 
dominant contribution to this process requires the presence of CP violating phases, and therefore 
we could expect to have a contribution in our model. As, due to gauge invariance, this decay 
depends on S'C/(2)i-breaking, it requires the presence of two Higgs vevs, either within 8lr insertions 
in the squark sector, or within gaugino-higgsino mixing on the neutralino and chargino sectors. 
At moderate tan/3, this disfavours the gluino contribution, since <5^ R insertions are, at most, 
proportional to the bottom mass. Of much more interest is the chargino contribution, as the 
(S^lr insertions are proportional to the top mass. Nonetheless, the required (S^ll and (6^i)rr 
insertions depend on powers of e, instead of e, which constrains greatly any chargino contribution 
dependent on the SU{3) structure. The chargino contribution will then be dominated by the 
flavour-changing of the CKM matrix, and thus cannot present any larger deviations than those 
predicted by MFV models. 



C. B physics 

Let us now discuss new physics effects on B s mixing. For this, it is convenient to adopt the 



following parametrization 



m-. 



M SM , M SUSY 

C Ba e ^ • (28) 

The B s -B s transition amplitude is divided into two parts, one arising from the SM loops and the 
other from the gluino-squark which are the largest SUSY contributions in our model: 

= (B s \H^m, = (B s \H^m. (29) 

Here, we focus on the phase <pB a rather than on Cb s , since the hadronic uncertainty in the latter 
is larger than the extra contributions that can be expected in this model. 

The current data of B s mixing phase is showing an interesting deviation from the SM prediction. 
A constrained fit, performed by HFAG, results in 

<Pb s = -0.36±g;i? or -1.171°;^. (30) 

Recall that non- vanishing (f>s s is an indication of new physics. The above fit is away from the SM 
at the level of 2.4 a, and the 90% CL range is: 

<t> Bs G [-0.61, -0.045] U [-1.48, -0.92]. (31) 
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It would be amusing to see whether or not our model could push this phase close to its best fit 
value. 

In order to estimate maximal size of 4>b s within this model, one should consider the ratio of 

Mfi 2 SY t° appearing in Eq. ([28]) . We follow Ref. [3] to express the supersymmetric amplitude 

in terms of mass insertions. The SM amplitude is available in Ref. 781 for instance. For this, we 

□ 

use Vts and Vtb from the full fit by the UTfit collaboration [79(. The result can be written in the 
form, 



^s S 12 SY _ 2if3 s 



0.71 2 0.062 2 



500 GeV\ 



2 



(32) 



where we have taken the ratio m~/m~ = 1, since we need only the order of magnitude of the above 
ratio. 

The first factor on the right hand side comes from the phase of Mf^ which is equal to — 2/3 s = 
-0.04. Now, from Eq. f[28|) . it is clear that a change in the B s phase, Arg [Mf™ + Mf^ 2 SY ] , would 
require \Mf^ 2 SY /Mf^\ ~ O(l) and a sizeable phase (/>|Y SY = Arg [Mf^/Mf^]. This fact has 
strong implications on the size and phases of the mass insertions appearing above. These mass 
insertions differ in different variations of our model, for instance, that largest Mis can be expected in 
RVV2. In this model, the mass insertions are (#,3) ll ~ £ — 0.05 and (S^rr ~ £2/°' 5 — 0.015\/ sec/3, 
as one can find in Eqs. (jlip . The size of the RR insertion depends on tan (3 and taking tan/3 = 10 
for example, we have (S^rr ~ 0.05. These values appear to be large enough to make a significant 
change in the B s phase from the second term of Eq. (|32p . Moreover, the product of these two 
insertions can have an O(l) phase. Notice that this is possible only in RVV2 where the leading 
terms in these mass insertions have phases. In Fig. [3l we show the different contours of </>s s . 

Indeed, in model RVV2 there is a region on the (mo,M 1/ / 2 ) plane where it i s possible to shift 



4>B a into the interval in Eq. (|31fl . at relatively high niQ and low M-y/2 [8O1 . [8l[. However, these 
regions of large SUSY contributions to (pB s are excluded by LFV and tend to have a too large 
contribution to ex as well. If we confine ourselves on a strip allowed by ex, the maximal value 
of 4>b s is around 10~ 4 . From the plots we can see that, in the case of RVV1, we are again in the 
same situation and for the strips allowed by ex we always have </>b 3 < 10~ 4 . Although these values 
could in principle vary due to the unknown O(l) coefficients in the soft terms, we can not expect 
variations in the order of magnitude. Thus, for low and moderate tan/3, our models are not able 
to provide a solution to the (j>B s anomaly, and, in this situation, we would expect this anomaly to 
disappear with the inclusion of further data 5 . If this is the case, our predicted deviation should 



The authors of [83j claim to be able to satisfy the "I>_g s anomaly and the ex tension in RVV2. We believe these 



25 




500 1000 1500 2000 

m 




500 1000 1500 2000 

m 



FIG. 3: Contours of cj> Bs for RVV1 (left) and RVV2 (right) in the m Q -M 1/2 plane for tan/3 = 10 and A = 
(top) and Aq = mo (bottom) . The color of an area indicates the size of 4>b s j which decreases as one moves in 
the order of dark to light blue. Each line between two adjacent areas is one of <pB, = 1CP 2 , 10 -3 , 10 -4 , 10~ 5 
(the contour corresponding to <pB, = —0.045, which would solve the discrepancy, is within the LFV excluded 
region in RVV2, and forbidden by direct bounds in RVV1). As in previous plots, green areas correspond to 
the currect LFV constraints and yellow area to the (g — 2) M favoured region. The area between the dashed 
black lines solve the ck tension, and the dark brown region show areas excluded by having a charged LSP 
or by LEP, excepting the Higgs mass bound, which is shown in thick dashed red lines. 



also be hard to observe at LHCb since it is much smaller than the precision of <f>B B attainable after 
5 years of run (10 fb" 1 ), which is estimated to be 5 x 10~ 3 [82]]. 

Similarly, we can repeat a parallel discussion on mixing with the following parameterization, 



CB d e 



2iAf3 



l dl2 



M SM 
lvl dl2 



(33) 



where each is the same as that in Eqs. (|29|) with the replacement s — > d. The supersymmetric 
to SM contribution ratio is given by 



^ll U 2 SY 



-2i/3 



3JLL 



(5 



d \2 
13 )RR 



(S d n)LL(5i 



13 )RR 



0.15 2 



0.013 2 



500 GeV 



(34) 



points correspond to a large value of tan (3 > 50. Moreover, we checked that the parameters can conspire to give 
$B a of O(0. 1), even for moderate values of tan/3. 
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One can get this result in the same way as for Eq. (|32p . 

In this case, each denominator in the square brackets is smaller than the corresponding one in 
Vtd/Vts\ 2 - Again, we use the CKM matrix elements from the full fit by the 



Eq. ((321) by the factor 



UTfit collaboration [7J]. In RVV2, we have (Sf 3 ) L L(St 3 )RR ~ ee 3 y£' 5 ~ 1.7 x 10- 5 v / sec^, and thus 
we could expect a sizeable effect on S^k if it were not for ex- However, AS^k is suppressed below 
1 x 10~ 3 on an sr strip in the same way as (fig^is. This is below the sensitivity of a super B factory 
to S^k whose estimate is around 5 x 10" 



601 ]. In RVV1, the new physics effect is even smaller, 
since the phases are suppressed. In the end, B^-Bd mixing is not very much affected in both 
models. This fact, a posteriori, justifies the way we determine the 0(1) coefficients in Yukawas: we 
tune the coefficients so that they reproduce the CKM matrix elements which were obtained under 
the assumption of no new physics. 

Let us come back to the alternative solutions of the ex puzzle that we mentioned in the previous 
subsection. If one is to blame the UT fit tension on S^k, one would need a change due to new 



physics of the amount AS^k ~ 0.07 [68j. However, this is two orders of magnitude bigger than 
what can be maximally expected in our models when e|^ SY is at the level of 18%, as we have 
seen above. In an alternative solution, e|F SY should be more suppressed, which in turn suppresses 
AS^k as well, thereby making it far less sufficient. The other possibility of invoking modification 
of AMb s I AMs d does not work for a similar reason. The authors of Ref. [831 ] need a new physics 
contribution to this ratio at the level of —22% in order to get an exact agreement. In our models 
however, the fractional change is about 0.5% at most, which is too small compared to what is 
needed. 

Finally, we can look to the expected values in our model for the decays b — > sfi + \i~ and b — ► svv 
that are similar to Kl — ► ifivv in the B sector. We find that the values of the branching ratios 
for these processes, in both RVV1 and RVV2 models for Aq = and Aq = mo, do not produce 
deviations from the SM values larger than the per cent level, in spite of the presence of sizeable 
flavour changing mass insertions. Therefore, these proceses are not interesting tests of new physics 
in our scenario. This is consistent with the work in 8J], where it was found that the deviation 



from the SM prediction is small once the constraints from b — > sj and B s — > /x + fx are taken into 
account. 
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D. Electric Dipole Moments 

The EDMs of fermions, such as the electron and the neutron, provide very stringent constraints 
on CP-violation in new physics. In the SM, the electron EDM from the Kobayashi-Maskawa 

The SM expectation for the neutron EDM is of 



phase is predicted to be < O(10 40 ) e cm 



O(10~ 31 ) e cm" 1 



, assuming a vanishing 6qcd- We can see that EDMs are highly suppressed 



in the SM, and thus they are excellent observables where to look for new physics. 



The electron EDM was studied in [4] within the context of RVV1. In these models CP is spon- 
taneously broken in the flavour sector. Therefore, the phases in the /i parameter and diagonal 
Af terms are very suppressed and can be neglected. In such a case, the imaginary parts required 



89)]. For the electron EDM, only 



for EDMs only appear from flavour-changing mass insertions 
neutralino-mediated diagrams contribute. The most important contribution to d e , when Aq = 0, 
comes from a bino-mediated diagram proportional to Qrn, [(813) ll {$33)1,11(831)111?], which is en- 
hanced by m T tan (5. 

In this work, we have updated our analysis on the electron EDM, taking into account the 
constraints on the $3 and x phases imposed by the fit of CKM matrix and the O(l) coefficients 
in Y e (see the Appendix for details). The results are shown for the RVV1 and RVV2 models in 
Figure HI An inspection of the slepton soft mass matrices indicates that, for RVV1, the leading 
phases in the product (<5f 3 ) l l{£ '33) L r{£ >fi) rr cancel and only subleading phases contribute, while 
RVV2 has non- vanishing leading phases. We then have the following predictions: 

(de)HWi ~ ^VbVt sm(2( X - AO) (35) 
(4)RW2 ~ ee 3 ^^ ' sin( X -2$) (36) 



901 ] and future 9l|, [92j d e experiments, for the 



In Figure H] we present the sensitivity of current 
two models. In this Figure, current constraints do not impose any restrictions on the parameter 
space of either RVV1 or RVV2. Nevertheless, electron EDM predictions are large enough to be 
probed at future EDM experiments. For relatively light SUSY masses we obtain d e ~ 10 -29 e cm -1 
and d e ~ 10~ 28 e cm" 1 , for RVV1 and RVV2, respectively. The latter predicts a value of d e about 
one order of magnitude larger than the former for any particular value of mo and M\/2 due to the 
larger e suppresion as seen in Eq. (|35p . This means that by reaching d e ~ 10 -29 e cm" 1 one could 
probe a much larger part of the evaluated parameter space, with tjiq < 1500 GeV, Mi/ 2 ~ 2000 
GeV. In particular, for RVV2, observation of SUSY at the LHC and solving the ex tension would 
force d e to be larger than 10 -29 e cm -1 . However, we have to take into account that these values 
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FIG. 4: Contours of |d e | = 1 x 1(T 28 c cm (dark red), |d e | = 5 x 1CT 29 c cm (light red) and \d e \ = 1 x lO" 29 e 
cm (brown) in the mo-M 1 / 2 plane for tan/3 =10 and Aq = (top), Aq = mo (bottom). We show predictions 
for RVV1 (left) and RVV2 (right). Current LFV bounds are also shown in green, and the (g — 2) M favoured 
region is shown hatched in yellow. The area between the dashed black lines solve the ex tension, and the 
dark brown region show areas excluded by having a charged LSP or by LEP, excepting the Higgs mass 
bound, which is shown in thick dashed red lines. 



will vary by factors O(l) because of the unknown O(l) coefficients to the different Mis. 

We now turn to the neutron EDM, d n . In this observable we have an additional difficulty, since 
its calculation as a function of partonic EDMs is not straightforward. Here we use two different 
approaches to this problem, the quark-parton 



931 ] and chiral quark 94j] models for d 



In the quark-parton model 93J, d n is expressed in terms of the up, down and strange quark 
EDMs, weighed by the fractional contribution A n of each quark to the spin of the neutron: 



dT = i! E (AX + A d n d d + A s n d t 



where rj E = 0.61 is a QCD correction factor 
A^ = 0.746 and A^ = -0.226. 
The chiral quark model 



951 ] and typical values for A' are A^ 



(37) 
-0.508, 



94| uses naive dimensional analysis to establish: 



d c ® 



(38) 
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where d q is the Wilson coefficient of the EDM operator at the hadronic scale. d q differs from d q , 
which is calculated at the SUSY-breaking scale, since the former receives contributions from the 
quark chromo-electric dipole moment, d^, and Weinberg's gluonic dimension-six operator, d G 



96, 



231 : 



96f|, and 



d q = r, E d q + V c ^-d c q + r, G ^d G (39) 

These contributions arise from the operator mixing ocurring in the running from the electroweak 
scale to the hadronic scale. In our estimates, we take the QCD correction factor rf = 3.4 
ignore d G , since it can be at most of the same order of d q and d q and we do not expect it to make 
big changes in our order of magnitude prediction. 

The (C)EDM of each quark has contributions from diagrams with gluinos d q , charginos d q and 
neutralinos d q . For gluinos and neutralinos, we can use the same arguments as in [4] to establish 
that the main part of each contribution will come from 3ro [(^3)^(^33) lr(^,i )rr] insertions. 
Chargino contributions come from both pure higgsino and wino-higgsino diagrams. The mixed 
wino-higgsino diagrams require no LR flip on the squark line. Contrary to the case with leptons, it 
is possible to have only one 5ll insertion, provided that the CKM matrix element Vij is not flavour- 
diagonal. For instance, the largest contributions to the down quark (C)EDM are proportional to 
Qto Vgd('5q„/)LLK* (i lrf . Pure higgsino diagrams require a LR flip on the squark line. One would 
expect this contribution to be strongly suppressed, as there are two Yukawa couplings on the 
vertices. However, it is possible to exchange one of the small Yukawa couplings for a yt or y^ 
coupling through a CKM off-diagonal term. As an example, the main contribution to the pure 
higgsino diagram for d^ is proportional to 3m [V^yt^^R^J^LL^yd]. Given the flavour 
structure of the SU (3) models, we find that the pure higgsino part shall be larger than the mixed 
wino-higgsino part, and will dominate the chargino contribution to d q c \ 



It was noticed in 



871 ] that, in addition to these contributions, it is important to consider 



beyond-leading-order (BLO) effects. These effects can be represented as new complex effective 



couplings [88], and although these are noticeable mainly for large values of tan /3, for EDMs they can 
also give significant contributions for moderate and small values of tan /3. In particular, diagrams 
with no imaginary part at LO can become complex from these effective vertices [891 ] . 



The most important BLO effect for the down quark EDM, for the values of tan (5 we are using, 
has been found to be the H contribution 



]. Although the BLO corrections also affect the gluino, 
neutralino, chargino and neutral Higgs loops, these contributions represent only small corrections 
and do not play an important role. Therefore, the only BLO effects we shall include will come 
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from the H± dia grams. 

The BLO effects can also enhance the coupling between the and fermions. At LO, the 
flavour suppression for an i/^-mediated dipole operator would be of order [ydV t * d y 2 Vu] ~ e 10 - 
The BLO vertex allows us to change the y d V^ d term for a (<5f 3 )_Rfj yb VA term, having then a total 
flavour suppression of order [(Sf 3 )RRybV t ly 2 Vtd] ~ e 6 - This is comparable to the gluino flavour 
suppression at LO. 

Even though the flavour suppression of the BLO contribution is of an order of magnitude 
comparable to the LO gluino contribution, the question that needs to be answered is how much 
does the extra loop suppression affect the diagram. This calculation has been done in for 
equal SUSY masses and moderate tan /3, giving: 

9m [V td (6 d 31 ) RR ] 



4 



a<i m 2 
9vr m 2 w 



2 

m 2 H± 



9m [OS&UlO&Jrh] 



9 

mr 



m 2 H± 



/ hSH (40) 



The first square bracket gives a factor of O(0. 1), which is compensated by the loop function 
f(m 2 /mjj±), which is of O(10). The third bracket is the flavour suppression ratio, which is of 
0(1). Thus, in this approximation, the importance of the contribution with respect to the 
gluino contribution depends mainly on the ratio between the squark and charged Higgs masses 
squared. Higgs bosons decouple differently from squarks and gauginos, as gluino and squark 
masses increase faster than the mass for increasing uiq and M]y 2 , the BLO contributions shall 
be more important for large values of these parameters. 

We shall include the contribution through the mass-insertion formulae of 



Such an 



approach contemplates three types of contributions. Two of them involve the CKM matrix, and 
are of the type 9m [V 3i [S^ll] and 9m [^((5^)^]. The third one is similar to that for the 
gluino loops: 9m [{^fs) LL^ii) rr] • The latter two shall be the ones that will allow us to avoid 
the y d suppression. The second term shall be particularly important for RVV1, where the gluino 
contribution has a cancellation between (5f 3 )LL and (S^rr- 

The SU(3) predictions for tan/3 = 10 and Aq = are shown in Figure [5] for the quark-parton 
and chiral quark EDMs, respectively. The EDMs receive a significant contribution from the 
loop in d d . Furthermore, dn P is also dominated by gluino and contributions to d s , while dn® 
is influenced by their contribution to d d , for both RVV1 and RVV2. The other contributions are 
smaller by at least half an order of magnitude. We can see that current bounds are always weaker 
than the LEP and LFV constraints and do not appear in the Figure. 

The observation of d n in the near-future experiments is not always compatible with the solution 
of the €k tension. In RVV1, which is mostly dominated by the loops, a d n of order 10 -29 
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FIG. 5: Quark-parton (top) and chiral quark (bottom) contours of \d n \ = 1 x 10~ 27 e cm (yellow) \d n \ = 
1 x 10 -28 e cm (orange) and \d n \ = 1 X 10~ 29 e cm (red) in the mo-M 1 / 2 plane for tan/3 = 10 and Ao = 0. 
We show predictions for RVV1 (left) and RVV2 (right). Current LFV bounds are also shown in green, and 
the (g — 2)^ favoured region is shown hatched in yellow. 



e cm" 1 is usually favoured, one order of magnitude under the reach of the next experiments 
In RVV2, both gluino and effects are comparable, the former dominating in the chiral-quark 
and the latter in the quark-parton models. In both situations it is possible to obtain an observable 
d n which is compatible with the solution of the ex tension. 

It is interesting to see that in some regions of the parameter space we have a cancellation. A 
careful analysis proves that close to these regions the phase of the (5f s )i,L term vanishes, and that 
can only be due to a cancellation between the initial term at Mqut and the contributions from 
the running. This causes a change of sign in the gluino contribution, such that it can interfere 
destructively with the part. We can avoid the cancellations by changing the sign of the O(l) 
terms at Mgut, m which case the interference is constructive. This, of course, will turn the 
cancellation into a small enhancement, but in any case, d n never exceeds its current bounds. 

Finally, Figure [6] shows the same information for Aq = thq. In this situation, the inclu- 



sion of the off-diagonal (Sfj)LR terms allow d n contributions of the type 9m (5fj)LL(<5ji)LR 



and 



(Sti)LR{6ti)RR 



These are always proportional to m&, although they do not receive a tan/3 
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FIG. 6: Quark-parton (top) and chiral quark (bottom) contours of \d n \ = 1 x 10~ 27 e cm (yellow) \d n \ = 
1 x 10 -28 e cm (orange) and \d n \ = 1 x 10~ 29 e cm (red) in the mQ-M 1 / 2 plane for tan/3 = 10 and Aq = too. 
We show predictions for RVV1 (left), RVV2 (right). Current LFV bounds are also shown in green, and the 
(g — 2)^ favoured region is shown hatched in yellow. 

enhancement. It was shown in [J] that these terms could be important, but they at most remained 
within the order of magnitude of the Sjm [(<5f 3 ) LL^is) Lii^fj) rr] terms. Notice that for d n the 
{$ij)LR insertions are of no interest: although RGE effects can make the off-diagonal terms as large 
as the (5fj)LR terms, they will usually enter the observables accompanied by a or {5^-) RR 

insertion, which are proportional to powers of e. 

An important contribution that appears when Aq ^ comes from flavour-diagonal subleading 
phases in (A^ii. Although the rotation to the SCKM basis removes the leading phases in the 
flavour-diagonal elements of the trilinear couplings, the different O(l) terms makes it impossible 
to remove the subleading phases in these terms. These phases are suppressed, but they still give 
important contributions, especially in RVV1, where the leading phases in the contributions from 
off-diagonal terms cancel. Such an effect is enhanced by the RGE evolution, and is particuarly 
relevant for the squark sector. 

In Figured] we can see that, barring the cancellations of the mass-insertion phases, both quark- 
parton and chiral quark EDMs are similar in magnitude. Although both RVV1 and RVV2 are still 
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not constrained by current bounds, the magnitude of the d n within the ex strip is now larger by an 
order of magnitude, enough to be observed in the near future. The main reason for the similarity 
in the order of magnitude between models lies in an enhancement of the gluino contribution, due 
to the (Ad) u subleading phase mentioned in the previous paragraph. For RVV2, this contribution 
mostly presents a correction to the off-diagonal contributions, but for RVV1 it becomes the main 
source for d n . 



IV. COMBINED ANALYSIS IN THE SUSY PARAMETER SPACE 

In this Section, we are going to present the results of a combined analysis of the observables 
studied in the previous sections. In particular, we will show the predictions of the model for LFV 
decays and EDMs in the regions of the parameter space where the SUSY contribution to ex can 
account for the possible tension between the measured value and the SM prediction, as discussed in 



641 ] . Such a requirement, and always up to possible variations of 0(1) coefficients, is fulfilled in a 
restricted portion of the parameter space and allows, as we will see, to do quite definite predictions 
for the other flavour observables. If we require in addition that the (g — 2)^ discrepancy between 
SM and data is explained by SUSY, we are restricted in a region of rather light SUSY masses, where 
most of the observables are expected to be close to the present experimental bounds. Given both 
the presence of unknown O(l) coefficients and the large theoretical uncertainties in the calculation 
which don't allow us to speak of a real failure of the SM, we cannot take what outlined above too 
seriously. Nevertheless, we think this can be a useful exercise in order to show how the interplay 
of various flavour observables can be used for testing this kind of flavour models. 

In order to understand the impact of the O(l) coefficients, in the following plots we shall set 
all of the 0(l)s in the soft mass matrices equal to unity. The 0(l)s of the trilinears shall be kept 
random, but fixed, since setting them to one aligns them with the Yukawas. At the end of this 
section we analyze the impact of varying the 0(1) coefficients, in order to understand how much 
they affect our correlations. 

As explained in Section IIIIB1 the SUSY contribution to ex depends on the sign of the 0(1) 
coefficient of the entry (6r R )i2, since a negative sign can induce a partial cancellation between the 
contributions proportional to (6^)12 and (6^)12(^1,1)12- Therefore, we have to take into account 
both the possible signs for (6j iR )i2- In Fig. [TJ we show the leptonic observables as functions of 
the lightest slepton mass for the model RVV2 with tan /3 = 10, Aq = 0, < mo < 2.5 TeV, 
< M1/2 < 1-5 TeV. The green band corresponds to the points of the parameter space for which 
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FIG. 7: Leptonic observables for model RVV2, tan /3 = 10, Aq = 0, as functions of the lightest slepton mass. 
See the text for details. 

(\ e K + e A" &y | ~~ l e ^ P |) < where we have checked that the correct sign has been obtained. For 



the SM contribution we took the central value computed in 63], | = 1.78 x 10" 3 , and a th is 
the corresponding theoretical error, cr th = 0.25 x 10~ 3 . The blue band represents the same for the 
case of negative (<5^ R )i2- The black and purple strips correspond to the more strict requirement 
+e^ sy | — l e /^ P |) < 3<r cxp , with cr exp being the experimental error, 0.012 x 10 -3 (in other words 



631]). The shaded regions 



for the strips we assume |e^| to take precisely to the central value in 
are ruled out by experiment, and the red lines indicate the future experimental sensitivity. For 
(g — 2)^, the area between the yellow lines solve the tension below 2a. 

The leptonic predictions in these "ej<--favoured" regions are quite interesting. We see that 
BR(/i — > ej) (top-left panel) gives at present practically no constraint, while the final MEG sensi- 
tivity (~ 10 -13 ), will test the model up to slepton masses around 0.5-1.2 TeV. Very interestingly, 
the sensitivity of a Super Flavour Factory reaching BR(r — ► ^7) ~ 10~ 9 can be rather similar 
(top-right panel). The reason for this is that the "e^-favoured" region selects rather low values of 
M1/2, where the future bounds of MEG and the Super Flavour Factory are comparable (see Fig. [T] 
in Section [III Aj) . Concerning eEDM (bottom-left panel), we see that reaching d e ~ 10 -29 e cm -1 
would test this case up to — 1.5 TeV, well beyond the reach of the LFV experiments. Finally, 
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FIG. 8: Leptonic observables for model RVV2, tan (3 = 30, Aq = 0, as functions of the lightest slepton mass. 

requiring that the SUSY contribution to (g — 2)^ (bottom-right panel) lowers the tension with 
the experiments below 2a, a rather light spectrum is selected, ^$ 250 — 500 GeV, so that all 
the other observables should be in the reach of running/future experiments, with, in particular, 
branching ratios of LFV decays being just below the present experimental limits. We checked that 
model RVV1 gives results for LFV decays and {g — 2)^ which are similar to the ones of RVV2, 
while d e is suppressed by approximately one order of magnitude: in this case, the future sensitivity 
on d e will test the parameter space up to slepton masses around 500-750 GeV. 

The predictions of Fig. [7] are rather stable also for oq = 1. In this case, both for RVV1 and 
RVV2, only BR(^ — ► cy) and d e gets slightly increased. For larger values of tan (3 (~ 30) we expect, 
both for RVV1 and RVV2, the future sensitivities for LFV decays and d e to completely test the 
parameter space in the range of parameters we considered, < rriQ < 2.5 TeV, < Myi < 1.5 
TeV. This is shown in Fig. [8] for RVV2. Notice that, also in this case, it is possible to account 
for the {g — 2)^ and the ex discrepancies at the same time, even if the leptonic observables are in 
general predicted to be very close to the present experimental bounds. 

In Fig. [9l we compare the discovery potential of the two most promising leptonic observables, 
fi — > e7 and the electron EDM. The correlation of BR(/^ — > ej) vs. \d e \ is plotted for both RVV1 
and RVV2, in the case tan/? = 10, ao = (left), tan/3 = 10, ao = 1 (center) and tan/3 = 30, ap = 
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FIG. 9: BR(/i — > &y) vs. \d e \ for different scenarios. See the text for details. 




FIG. 10: Lcptonic observables for model RVV2, tan/? = 10, Ao = 0, as functions of the lightest slepton 
mass, for a random variation of the 0(1) coefficients (red dots). See the text for details. 

(right). As before we studied the mass range: < rriQ < 2.5 TeV, < M\i2 < 1.5 TeV. In the 
figures, only the "e^-favoured" region with positive (<5^r)i2 has been plotted. The horizontal line 
corresponds to the final sensitivity of MEG, the vertical line to the sensitivity on \d e \ of the running 
Yale-PdO experiment. We see that, for RVV1, \i — > e^f should be able to constrain the parameter 
space more strongly than eEDM, while for RVV2 it is \d e \ the most sensitive observable (except for 
the large tan/3 case). These features could be useful in the future, in order to discriminate among 
different models and, more in general, shed light on the structure of mixings and phases in the 
slepton mass matrix. 

Finally, let us briefly comment about the impact of the unknown 0(1) coefficients on the analysis 
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FIG. 11: Same as Fig. [51 but with a random variation of the 0(1) coefficients (red dots). See the text for 
details. 

outlined above. In Fig. [TU1 we plot the same observables of Fig. [7J (RVV2, tan (3 = 10, a = 0) 
in red, performing a random variation of all the 0(1) coefficients between 0.5 and 2 (in absolute 
value). For comparison, we superimpose to the scatter plot the green and blue bands of Fig. [JJ We 
do the same for Fig. [9] in Fig. [TTJ only for RVV2, in yellow, superimposing the green bands. 

From these Figures we see that, although the lines are broadened by the variation of 0(l)s, 
as expected, the correlations are roughly mantained. In some cases we can see variations of even 
several orders of magnitude for the predictions. This is due to two different sources. On one hand, 
for fixed values of mo and M\i2 we can have a variation of a factor 4 or 1/4 in the product of two 
sleptonic 0(1) coefficients. On the other hand, if we allow the variation of O(l) coefficients in ex, 
the lines in the (mo, M1/2) plane of Figure [2] become broad bands, increasing in turn the width of 
the bands of the sleptonic observables. Still, our predictions remain stable enough, especially for 
some observables such as BR(// — ► e"f), to conclude that the main qualitative features discussed 
above are not affected too much by the unknown 0(1) coefficients. 



V. CONCLUSIONS 



In this work, we have studied the phenomenology of a MSSM with a SU(3) flavour symmetry 
and spontaneous CP violation. We have shown that in this framework it is possible to fit the 
observed fermionic masses and mixings simultaneously solving the so-called SUSY flavour and CP 
problems. As a proof of existence, we have analyzed an explicit example based in the model of Ross, 
Velasco and Vives [3] that successfully overcomes the present FCNC and CP-violation constraints. 
At the same time, this model predicts a non-trivial flavour structure in the SUSY soft-breaking 
terms that will show up in near- future CP-violation, LFV and hadronic FCNC experiments. We 
have analyzed a model that includes the minimal set of contributions to the soft-breaking terms 
and we have presented two possible variations with additional contributions to soft-breaking terms. 
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In these models, we naturally expect to be able to observe the \i — > cy decay at MEG for 
sfermion masses within LHC reach. Similarly, it is also possible to measure the r — ► fij branching 
ratio at the future Super Flavour Factory. In the hadronic sector the main effect of the model is a 
sizeable contribution to ex that could explain the recently observed discrepancy between the SM 



and the measured value taking the latest lattice values for the Bk factor [63|. The electric dipole 
moment of the electron is also probably to be within reach of the present experiments while the 
neutron EDM will be just beyond the planned sensitivity of the experiments. It is unfortunate 
that, at low tan/? and without any conspiration of 0(1) parameters, the evaluated versions of this 
model cannot provide the required contribution to (f>B s in order to solve the discrepancy reported 
in [75]. This means that if the future data confirms the discrepancy, non-minimal versions of the 
model shall be required. 

We have made a combined analysis of all the observables if we require that the ex discrepancy is 
solved by the SUSY contributions of the model. In this case we are able to relate the values of the 
different observables and again \i — > &y, r — > /U7 and d e can be measured in the future experiments 
if SUSY is to be found at LHC. 



Acknowledgments 

We would like to thank G. G. Ross, I. de Medeiros Varzielas and M. Passera for useful discus- 
sions. L. C. wishes to thank the University of Valencia for the kind hospitality and partial support 
during various stages of this work. J. J. P. thanks the universities of Padova and Wiirzburg for 
hospitality during his visits to complete this work. A. M., J.-h. P. and O. V. would like to thank the 
Yukawa Institute, Kyoto and the research project YKIS for providing the opportunity to develop 
an important part of this work in a stimulating and active environment. We acknowledge support 
from the MEC-INFN agreements FPA2008-04065-E and INFN2008-016. The work of O. V. and 
J. J. P. was supported by the Spanish MICINN and FEDER (EC) Grant No. FPA2008-02878. 
O. V. was supported in part by European program MRTN-CT-2006-035482 "Flavianet" and the 
Generalitat Valenciana under the grant PROMETEO/2008/004. A. M. acknowledges partial sup- 
port from the CARIPARO Project of Excellence (LHCosmology). W. P. is partially supported by 
the German Ministry of Education and Research (BMBF) under contract 05HT6WWA. J.-h. P. ac- 
knowledges Research Grants funded jointly by the Italian Ministero dell'Istruzione, dell'Universita 
e della Ricerca (MIUR), by the University of Padova and by INFN within the Astroparticle Physics 
Project. A. M. and W. P. are partially supported by MRTN-CT-2006-035505 (Hep-Tools). A. M. 



39 



and J.-h. P. acknowledge partial support from the INFN FA51 Grant and from the RTN European 
Contracts MRTN-CT-2006-035863 (UniverseNet). 



APPENDIX A: CKM FIT 

The SU(3) model generates the following leading structure for the Yukawa matrices, in terms 
of e and e: 



Y d 







i xf 3 e 3 e iSd+/33 x% 3 S d e 2 e if5z 



x 13 e e 
x d 23 Z d e 2 e^ 
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where 5f = 2a f + 0% + f3' 2 ■ It is necessary to fix the values for the xfj parameters and the flavon 
phases. In the quark sector, this means that we need to reproduce the quark masses and CKM 
matrix and, in fact, we will see that the CKM matrix is the main source of constraints on xfj and 

After diagonalizing the Yukawas, keeping terms up to order e 3 , e 2 and ee 2 , and rephasing fields 
such that V u d, V us , V cs , Vd, and Vtb are real, we get the following CKM matrix: 



Vckm 
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where A a are complex corrections due to subleading terms, with effective phase to a . These sub- 
leading terms play a very important role in the CKM phase 5ckm which consists on a combination 
of subleading phases, the largest of these being u) u \,: 



followed by uj us : 
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Other A Q coefficients are A u d, A c fe ~ l+0(e 2 ) and do not play an important role in the fit. However, 
it is important to point out that the structure of these subleading terms is highly model dependent, 
and that the values we get for them shall be valid only in this model. Other SU(3) models can 
have a completely different subleading structure. 



The phenomenological fit of 36j, based on structures similar to those in Eq. (lAlft . was used in 
previous works 3, 4] to fix the 0(1) parameters of the Yukawa matrices. In this work we make 
a new fit specifically for the SU(3) model, taking into account also flavon phases. Although the 
variation of the latter make very small changes in the CKM elements, such variations could make 
these elements exceed the 3a bounds in 3?|]. 



The values found in [3jJ for the CKM matrix are connected to a large number of flavour 
observables. The analysis assumes that the SM is the only source of flavour and CP violation. 
Nonetheless, the strong consistency in the SM between all flavour observables, in particular those 
participating in the construction of the unitarity triangle, suggests that new physics effects are 
small, and that the main source of the registered flavour and CP violation still lies in the CKM 
matrix. Thus, we shall adjust the O(l) parameters in Y u and Yd such that the CKM matrix of [stJ 
is reproduced, assuming that SUSY contributions do not affect the fit significantly. 

To make a fit of the 11-dimensional parameter space of Eq. (|A1|) . the Powell minimization 
method was used. The function to minimize was defined as a x 2 on the quark masses and on the 
four CKM parameters in the Wolfenstein parametrization: A, A, p and fj, as in j^. We required 
the 0(1) parameters not to be larger than 2, and not smaller than 0.4. We found that the following 
values give a very good fit on the masses and mixings: 
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361 ]. with 



with P' 2 , /?3 and x non-uniquely determined. The results for xfj are similar to those in 
a larger discrepancy in x 2 ^ • The value of (a u — ad) is mainly fixed by the CKM phase and | V us \ , 
as reported in Qj. To leading order, we can extract a la uncertainty afj on the xfj constants 
from the errors on the CKM parameters We roughly expect af 2 ~ 7 x 10 -3 , af 3 ~ 5 x 10~ 2 
and a"23 ~ 4 x 10~ 2 . As the xfj parameters participate only in subleading terms, such a rough 
estimation is not possible. However, by fixing these constants at the above values, we can estimate 
the uncertainty on a u — ad, being ar au _ ad \ ~ 0.04. 
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FIG. 12: Fit to Quark Masses and CKM matrix at one (red), two (blue) and three (green) sigma. There is 
a strong correlation between /? 3 and \i but (3' 2 does not affect CKM elements. 





m u 


m c 


m t 


m d 


m s 


m b 


SU(3) (GeV) 
Reference (GeV) 


2.9 x 1CT 3 
(1.4 ±0.5) x 1CT 3 


0.57 
0.63 ± .08 


172.0 
170.3 ±2.4 


4.1 x 10~ 3 
(3 ±1.2) x 10~ 3 


71 x 10~ 3 
(56 ±16) x 10~ 3 


2.85 
2.89 ±0.11 



TABLE VII: SU(3) predictions for quark masses at n%z and corresponding values obtained by running to 
the same scale 1991. 



After fixing the determined 0(l)s at these values, and taking ad = 0, we make a grid-based 
X 2 analysis on the remaining three phases. We find out that the CKM matrix is completely 
independent of (3' 2 , as expected from Eq. (|A3p . The result of the scan for /3 3 and % is shown in 
Figure [F2l These two phases enter through subleading terms in the CKM matrix, and are severly 
constrained due to the precise determination of \ V US \, |V^,| and the diagonal elements. We are thus 
forced to place (03, x) around (—20°, 20°), or on any of four degenerate spots obtained by adding 
180° to each. For definiteness, in the following we set (/3 3 ,x) = (-20°, 20°), and (3' 2 = 0°. 

Using these values, we compare the results of the fit with the measured values for the CKM 
matrix and quark masses at mz [99(] . The results are shown in Tables IVIII and IVIII1 We can see 
that the fit is very good and almost all of the SM parameters lie within their la values. 

Regarding the lepton sector, we shall unify the structure of quark and lepton flavoured matrices, 
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A 


A 


P 


V 


SU(3) 
Reference 


0.2250 

n 99^7+0.0009 
u.zzo i -o.ooio 


0.805 


0.156 

n ioc+0.031 


0.327 

U.«5*y_ .017 



TABLE VIII: SU(3) predictions for CKM elements and corresponding measured values, taken from 37 1. 





m e 




m T 


SU{3) (GeV) 
Reference (GeV) 


4.867808 x 10~ 4 
4.866613(36) x 10~ 4 


104.47907 x 10~ 3 
102.72899(44) x 10~ 3 


1.77634 

L.<10tO_ .00026 



TABLE IX: SU(3) predictions for lepton masses at mz and corresponding values obtained by running to 
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with and S e = and T, u = 0. Notice that in Y u the zeroes are actually terms of order higher 
than e 3 . We use the same 0(l)s in Y u for Y v and we find that the 0(l)s that fit best the charged 
lepton masses are: 



/ xf 2 \ / 1.31 \ 



x\ 
xS 
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(A9) 



V 2 -° J 



13 

\ x h j 

which give the values in Table IIXI Note that lepton masses are known very precisely and thus 
subleading contributions will play an important role in the fit. For this reason, we did not try to 
fit exactly the reference values in Table IIXI 

For more details on how the SU(3) model accomodates neutrino mixing, we refer the interested 
reader to 



jie 

a. 
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APPENDIX B: O(l) COEFFICIENTS IN THE SCKM BASIS AT THE GUT SCALE 

In general, flavour symmetry models contain unknown O(l) coefficients at the flavour scale. 
In CPV and FCNC studies we are interested in the leading real and imaginary contributions. 
However, we find that in some particular observables, leading phases might cancel. In such a 
situation it is important to know the structure of subleading complex terms in order to make 
verifiable predictions. For this purpose, we present here the structure of the O(l) coefficients for 
the leading terms after the SCKM rotation. 

In the following, we shall define x$ = (x{ 2 x 23 ~ x { 3 )' with f = u,d, e. 



1. RVV1 



M 



Q 



Our soft matrices in the flavour basis have the following structure at the GUT scale: 

/ 1 + Q 1 e 2 y t N 

l + Q 2 e 2 Q 3 e 2 e^ 

\ Q 3 e 2 e-^ 1 + Q 4 y t j 



ml, 



with exactly the same structure, but different 0(1) constants Ui and Lj, for M? c and M~, respec- 
tively. For M 2 - and M~ c we have an analogous structure, with e replaced by e, yt replaced by y b , 

d R r 

and different O(l) constants Di and E^. This shall be true for all variations of the model, so in 

following subsections we shall only show M? without repeating these specifications. 

w 

Rotation into the SCKM basis gives Eq ([9]), with the following O(l) coefficients: 
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It 

M% (3, 2) = D 3 - D A xf 3 E d y b e" 2 ^^) (B2c) 
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M?(2,l) = (Q 2 -Qiyt)(xi 2 /Z d ) (B3a) 
M?(3,l) = -Q A xj + ^9l {x f 2 /j: d)e Hx-M (B3 b) 

M? (3, 2) = Q 4 4, E d - e nx-M (B3c) 



M|c(2,l) = {E 2 -E 1 y h ){x\ 2 /T >e ) (B4a) 
M|(3,l) = ^(x^/Se)-^^^^- 2 ^-^) (B4b) 
Ml (3, 2) = E 3 - E A x e 23 S e y 6 e ~ 2 ^~^ (B4c) 
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RVV2 



The Kahler potential in RVV2 has got the minimal terms with the addition of the #3^23 flavon 
vevs. This modifies the (2, 3) and (3, 2) elements of the soft mass matrices, giving them the 
following structure at the GUT scale: 

f l + Q l e 2 y t \ 

M|= 1 + Q 2 e 2 Q 5 ey^e^ 

\ Q 5 ey^e-^ 1 + Q 4 y t ) 

For simplicity, we have omitted the minimal X 3 (X = U, D, Q, E, L) 0(l)s, replacing them by 
the corresponding X§ that accompany the #3 #23 flavon vev. The O(l) coefficients are: 
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M% (2, 1) = {D 2 - Dt y b )(xi 2 /X d ) (B7a) 

M?(3,l) = D 5 (x d 12 /Z d ) - eD 4 x d s y° b 5 e- 1 ^ 3 -^ (B7b) 
"a 
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Mi (3,1) = E 5 (x e 12 /J: e )-eE A x e s y b 5 e- i ^-^ (B9b) 
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3. RVV Model 3 



The addition of an effective term with the antisymmetric (#3# 2 3)#3 flavon vevs gives the soft 
matrices the following structure, in the flavour basis, at the GUT scale: 
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The additional term in the (1,3) and (3, 1) sectors of the mass matrices induce large terms in 
the (1,2) and (2, 1) sectors of the mass matrices, and have a moderate influence on the (2,3) and 



(3,2) sectors. The structure of the 0(1) coefficients follow: 

Ml =(2,1) = (U 2 -U 1 y t )(x^/Z u ) + U 6 x% 3 Z u y t e- i ^ +5 ^ (Blla) 

N. 

M|c(3,l) = U 6 (Bllb) 

M4(3, 2) = U 3 - U 4 a&E u yt e 2i ^ + U 6 y t {x u n /T, u ) e^-M (Bile) 

M% (2, 1) = (D 2 - Di y b )(xj 2 /E d ) + D 6 x d 23 E dI/6e *(2x-fc-*<0 (B12a) 

Ft 

M? (3, 1) = As (B12b) 

a R 

M% (3, 2) = D 3 - D 4 x d 23 S d y b e -*<X-to) + D 6 y b (x? 2 /E d ) e l ^~^ (B12c) 

M? (2, 1) = Q 6 x^ 3 S d + £ (Q2 ~ Qiyf) (^ 2 /E d ) e *H*+A+^ (B l3a) 
y £ yt 

M?(3,l) = Q 6 (B13b) 

M|(3,2) =Q 4 x^ 3 S d -| Q 6 (xf 2 /S d )e- W " + —e"** e* 2 *"^ (B13c) 

Mi (2,1) = (E 2 - Ay b )(x? 2 /S e )+ J B 6 xl 3 S e y 6 e i ( 2x - A -« (B14a) 

-TL 

M| c (3, 1) = E e (B14b) 

M|c (3, 2) = £ 3 - £ 4 *I 3 S e y b e - 2i (x-^) + ^ yb ( x ^ 2 /s e ) e ^~« (B14c) 



M|(2,l) 
Mf(3,l) 
M|(3,2) 



Lq x 23 S e + 



L4 x| 3 S e 



e (L 2 - Liy t ) 



lit 



(xf 2 /E e ) e 



i(-2 X +/3 3 +5 (i ) 



L 6 (xf 2 /S e )e" 



. £L 3 _ 
+ - — e 



3 i(2x-/3 3 ) 



(B15a) 
(B15b) 
(B15c) 
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A- Terms 



For all of the evaluated models, the A- Terms will have the same structure the Yukawas have, 
but with different O(l) constants. We parametrize them in the following way: 

( Af xf 2 e 3 e i5d A d 2 xf 3 e 3 e i ^ 3+5 ^ ^ 

A d = A d x( 2 e 3 e i5 J A$Z d e 2 Ajx^Zje 2 e^ 3 A y b (B16) 

y A^x^e 3 ^^ 5 ^ Af x$ 3 S d e 2 e ift Afe 2 ^ j 

and similarly for A u and A e . 

The A- Terms are rotated into the SCKM basis using the same unitary matrices that diagonalize 
the Yukawas. However, as the 0(1) constants in each entry are different from the ones in the 
Yukawas, the A- Terms are not diagonalized, and mantain the same structure. For squarks, the 
remaining off-diagonal terms are then affected by the rephasings, and we obtain: 

/ n A fu x u^ £ 3 e -iw' A fu £ 3 e i(2 X -u/) \ 

A 3 S u e 2 



Af x% 3 S u e 2 e 2i * 



\ A' 2 U £ 3 e i(^'+2/3 3 -2x) A 'u ^ ^ £ 2 g2i(/3 3 -x) 



/ 



A n 



A' d 



Id „d =3 „— iuus 
12 E e 



A'-fx 



A3* 



4" 
^5 



x 13 6 e 



(B17) 



A£* xf 3 e 3 e *(^-+2A-2 X ) ^ ^ Sd ^2 e 2i(/3 3 -x) 



^■4^ x 23 ^ 2 



A y b (B18) 



/i/e e =2 
J±2 x 12 fc 





Ale ™e ? 3 

^1 x 12 fc ^3 1 

\^ A' 2 e x\ 3 e 3 e^ 3 "*) x\\ S e e 2 e^ 3 "*) 



g3gi(/3 3 -x) \ 



A| S R £ 2 



Ale e 
^2 x 13 



Af x| 3 S e e 2 e*C^-x) 
^5 



A y& (B19) 



with: 



a/ = (4 - 4) 

4'/ _ I if _ /\f\ _ ( _ Af\(r r f ryf InfS \ 

^2 ~~ 1^2 A b> ^5 A x 12 x 23/ X 13^ 



.4 



-A 



(B20a) 
(B20b) 
(B20c) 



with no subleading phases generated by the SCKM rotation. Notice that if the initial A{ = 1, the 
A-terms are aligned with the Yukawas. Thus, the new A'/ go to zero, and the A-terms are diagonal 
at the SCKM basis. 
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APPENDIX C: O(l) CONSTANTS DUE TO RGES 

Using the leading log approximation, we can estimate the effect of the running on the flavoured 
matrices. For this we define the negative parameter A, such that: 

> O(l) 



1 1 16vr 



i MsuSY 

log 



M G UT 

The rotation to the SCKM basis must be done at Mw- However, since the off-diagonal RGE 
contributions to the Yukawa matrices are negligible, we can perform the rotation at the GUT 
scale, and then apply the corrections from the running. 

For M^, one must take into account that only the (3, 3) term of Y v gives a significant contribution 
to the running. This parameter only participates until the heavy (z^)3 decouples, at the scale M3. 
To take into account this effect, we introduce the parameter: 

A 1 1 1 M 3 



A 167r Mqut 

With the exception of the A- Terms, the general result for all models is that the structure in 
terms of e and e remains unchanged. Even though the generated parameters are not necessarily of 
O(l), the A parameter provides an additional suppression, which approximately keeps the matrix 
structure very similar to the one at the GUT scale. 

These induced coefficients can have two parts: one independent of the flavour structure of the 
soft masses (the MFV contribution, generated by the misalignment of the Yukawa matrices), and 
one dependent on the flavour structure. In the following equations, those parts independent of the 
soft mass flavour structure can be distinguished because they do not vanish when the Xi terms 
{X = U, D, Q, E, L) are set to zero and when the A, terms are set to one (this also sets all primed 
variables to zero). They are only present in the M? M~ and Ai matrices. 

In the following, we shall denote Mj~ (z, j), at the scale /j,, as X^, with X = U, D, Q, E, L. For all 
soft mass terms, we shall factorize the e and e factors that appear in the mass matrices presented 
in Eq. Q, as well as the leading phase. Since we shall also factorize tbq, we define ao = Ao/mo. 

1. RVV Model 1 



U£ w = U£ UT (Cla) 
Ui w = Ug UT + 2{u£ UT -2atA?A% x ^}ylA (Clb) 
Ug w = Ug UT + 2{ug UT + 2alAtAl X ^ u e^}y^ (Clc) 
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D™ = Dg UT (C2a) 
Dg W = Dg UT + 2 [Dg UT - 2a 2 A' d A d x d 3 e -*<x-AO } y 2 A (C2b) 
Dg W = Dg UT + 2 + 2a 2 ^ Af x d 23 S d e -*<X-A) } y 2 A (C2c) 



Qfi W = Qg UT -2^(3 + a 2 (A u 5 ) 2 )x d x d 23 Z d y 2 A 

{2(Q 4 ft + U 4 y t )x d x d 23 y 2 t + 2a 2 (A ,d A' d x d 13 x d 23 + ^Afx?^ 2 } S d A (C3a) 
Ql W = QfF T ^-2(3 + Gg(^) 2 )x^ t A + 2^(3 + a2 J 4^)( 3 ;f 2 /S d K3E u 2/ t e 2 ^A 
- {QsT T (% 2 + </?) + (Q4 + 2C4)4y? + 2ag 4? xf 3 | } A 
+ £ ^(Q4 + 2U 4 )(x 12 /Z d )x 2 l 3 Z u y 2 e 2i XA (C3b) 
Qf 2 W = Q^ T + 2(3 + G 2 (^) 2 )43S,y t A-2^(3 + a 2 ^K3S uyi e 2 ^A 

+ I^V^ 2 + % 2 ) + (04 + 2C/ 4 ) f 1 - g^e 2 ^) y 2 - 2a 2 H 4 4 ^j x d 3 ^C3c) 

1^23^ V e ^23^0! / Mi J 



^ = ^T T (C4a) 
£ 3 T = ^ T + 2{^3T T -2a 2 ^^ 5 xf 3 e- 2 ^-*)}, 2 A (C4b) 
Eg W = Eg UT + 2 {Eg UT + 2a 2 A'i A% x e 23 Z e e~ 2 ^-^} y 2 b A (C4c) 



T EW _. r GUT A — T e -r e y I, 2 A A 
^21 — ^21 ~~ 4 ^2 X <5 X 23 2/t ^ 

{2(L 4 y t )x% x\ 3 y 2 A, + 2a 2 (A' 2 e A'ixi 3 x e 23 + ^ e ^x? 2 )y 2 } S e A (C5a) 
■ EW _ nGUT 



L 3 \ = Q^ UI - 4xg y t A v A 



Qaf 1, (% 2 + yt A„) + L 4 ^ 2 A v + 2a 2 A'i A% x\ 3 V j- } A (C5b) 



Lf 2 W = Qs2 UT + 4x e 23 S e yt A, A 

r o GUT v 2 

+ ^§V(y 2 + Vt A,) + W A, - 2a 2 4f^ } x e 23 y e A (C5c) 

i ^23^6 yt 
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2. RVV Model 2 



ttEW 
u 21 



ttGUT 
u 2l 



u, 



EW 
31 



Ug UT + 2 { Ug UT - 2al jig A? A- xfi, e^+ft) } y 2 A 



U. 



EW 
32 



0.5 4 ^5 ^23 w e 



(C6a) 
(C6b) 

(C6c) 



D 2 T = ^ 2 T T (C7a) 
= ^ T +{2^ 3 T T -2a^A' 2 ^^f 3 e-^-^} y6 2 A (C7b) 

= ^T T + 2{^ 3 T T + 2a^<^^3^e-^^)} y 2 A (C7c) 



Qf^ = Qg UT + 2^(3 + at(Atf)xlx d 23 X d y?A 



e 4 



+ 

=r2 



{2(Q 4 ift + 1/4 4s Vt + 2a 2 (A 2 d ^xf 3 x^ + A[ d Aixi 2 )y 2 b } S rf A 
Q 5 % a5 (xf S-^^) + 4 2 4 3 e ^A-ft)) y 2 A (C8a) 



e 
e 

Q3T = Q 3 7 r + 2^(3 + a 2 (^) 2 )4^ e 4( - 2x+/33+ ^ ) A 

Q 3 T T (% 2 + Vt) - - ( (Q* + 2^4)4%°- 5 % 2 + ^^-^ e^+ft+^l A 



(C8b) 



£ Vt 



+ {Qg UT (y? + Vl) - j ((Q4 + JH/^V + 2alA'tAl^j X d 23 E d e^^ } A 

(C8c) 



'21 


pGUT 








'31 


- ^31 


+ { 


2^T r 


2fl ,.0.5 

yb 


i£W 
'32 


pGUT 
- ^32 


+ 2 




+ 2a "OB" 

yb 



(C9a) 



ytA (C9b) 
^ e AI S e e -^x-A-^) | y 2 A (C9c) 
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+^ {2(L 4 yt )x% x%% y 2 t A v + 2a 2 (A'£ A% x\ 3 x\ 3 + Af A e 3 x e 12 )y 2 } S e A 

Jl L5 y0.5 e i { -2 X+M > 2 ) + j^-fr-f^ y 2 A ^ A (C1Qa) 

Lff = L^ + ^^yfe^+^+^A.A 

- {Z 3 T T (Z/| A, + y 2 ) - f (l^V?" 5 !/? A, + 2a 2 ^z? 3 J^) e «-**A+fc) | A 

(ClOb) 

£ Vt 

+ \l^{uI a v + y 2 ) - 1 [l^-SI a, + ^Iaui^e) x^-^+^a 

(ClOc) 

3. RVV Model 3 



U£ w = U§r (Clla) 
U£ W = u£ UT + 2U£ UT y 2 A (Cllb) 
^ = C/ 3 G 2 ^ + 2{^ + 2a 2 ^S ue 2 ^} y 2A (Cllc) 



£>Jf = Dg UT (C12a) 
= Dg UT + 2D^ T ylA (C12b) 
D 3 T = ^sT T + 2 {D 3 T T + 2a 2 , AfAl S d e -*<*-*> } ^ A (C12c) 



QT = Q^ T + (l + ^ ^e" 3 *) (Qe^rfA (C13a) 

QfT = + Qe(y 2 + yl)A (Ci3b) 

Of = Q^ T + 2(3 + a 2 (^) 2 )4 3 S,y 4 A-2|J(3 + a2^)43S M y t e 2 ^A 

+ { + Vl) + + ^) (l - f^e 3 *) , 2 - 2a 2 Af A'f V j ) x^ d A 

(C13c) 
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Eg w = Eg UT (C14a) 

E EW = E gUT + 2E GUT y 2 A (CMb) 

Eg W = Eg UT + 2 {Eg UT + 2al A^A% e" 2 ^)} y 2 A (C14c) 



Lg w = Lg UT + (L 6 x e 23 S e )y 2 A, A (C15a) 
Lg W = Lg UT + L 6 ( y 2 A, + y 2 )A (C15b) 

( tGUT 2^, 
+ "tV ^ + % 2 A ^ + W A - - 2fl A 5 ^ - ^ 3 E e A (C15C) 

4. A- Terms 

To write the RGE contribution for the A- Terms, we shall first denote: 

4 = ~3i ~ 3g 2 2 - y gl + 3y 2 + y 2 (C16a) 
«d = ^gfM, + 3<? 2 M 2 + y 5 |M 3 + (3y 2 b Af + y 2 A e 5 )A) (C16b) 



-I Q "1 /? 

k£ = — 5 ?M! + 3<? 2 M 2 + y 5 2 A/ 3 + 3y?A%Ao (C16d) 

«e = - 3 5 1 + 3y 2 + y 2 (C16e) 

5 

k\ = ^-glM 1+ 3glM 2 + {ZylAi + y 2 T At)A (C16f) 

As the A- Terms are not hermitian matrices, the RGE evolution is different for each off-diagonal 
term. For each element, we shall factorize in Eq. (|Bf 7|) the leading order of magnitude in terms of 
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e or e, the phase, and the global term A^i/i. Then, after the running, we get: 

A d (l,l) = Oie 4 ) (C17a 

A d (l,2) = (l + 4/\)A'fx d l2 (C17b 

4,(1, 3) = (l + (k£ + A*™ 1 ) a) A' d xf 3 (C17c 

A d (2,l) = A d (l,2) (C17d 

A d (2, 2) = (l + (k$ + 2^£j a) A* S d (C17e 

A d (2, 3) = (l + (4 + A^ m ) a) Af x d 23 S d (C17f 

Ai(3,l) = (l + (^ + 5y 2 + y 2 )A)A' 2 d x? 3 (C17g 

A d (3,2) = (l + (Kj + 5y 2 b +y^)A)A' d x d 23 ^ d (C17h 

3) = (i + (% 2 + + 4 + a) ^ (cm 

(ci7j; 



with: 



.4 



d,run 
13 



,4 



d,run 
23 



.4 



13 



A 



23 



4*(M 
4,(1,2 
4.(1,3 

4*(2,1 
A u (2,2 

A u (2,3 
4,(3,1 
A u (3,2 
A u (3, 3 



0(e 4 ) 

(1 + («J* + A™)A) A^ 3 
4.(1,2) 

1 + 04 + A™) A) Afx u 23 E u 
l + (^ + 5y 2 + y 6 2 )A) A^ 3 
l + (^ + 5y 2 + y 2 )A) Af^ 3 S u 



A I Al 



' O A u 4- A d \ T d / 9 A u 4- A d \ f T d r r u V \ 

2 / 13 V 2 / ' X-\nZ-i r , 



a 2 o ^5 2 , ff2 { 2A% + Af \ ( x 23 T, u \ 



%t + Is 



2Ai + A^ 



4/U 
A 2 



x 13 t(w'-w„ a ) _ fl 



13 



e 



A'u 



13 u 



j4 4 / fc y j4 4 / \x 23 Zj u 



(C18a) 
(C18b) 
(C18c) 
(C18d) 
(C18e) 
(C18f) 
(C18g) 
(C18h) 
(C18i) 

(C19a) 

(C19b) 

e- 2 ^C19c) 

(C19d) 
(C19e) 
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We see that the running can generate subleading phases, which are contained within the A('^ un 
terms. In particular, it is noticeable that these phases do not vanish when aligning the A- Terms 
with the Yukawas (setting all A'/ — > 0), which means they belong to a MFV contribution. 

It is also important to remark the fact that A\o Un and A^ un contain terms of order (e/e) 3 
and (e/e) 2 , which are enhancement factors. Again, these are due to MFV contributions. Thus, the 
RGE evolution has the potential of changing the structure of A u noticeably. 

We do not list the A e matrices, since their structure is identical to A^, with the replacements 
Vt Vt i^d K e ) x f x i ) M —* A\ and A^ — > 0. One must also consider S e = 31^ and 
E„ = 0. 
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